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ABSTRACT 

Atmospheric parameters and oxygen abundances of 825 nearby FGK stars are derived using high- 
quality spectra and a non-LTE analysis of the 777 nm Oi triplet lines. We assign a kinematic prob- 
ability for the stars to be thin-disk (Pi), thick-disk (P2), and halo (P3) members. We confirm pre- 
vious findings of enhanced [O/Fe] in thick-disk (P2 > 0.5) relative to thin-disk (Pi > 0.5) stars with 
[Fc/H] < —0.2, as well as a "knee" that connects the mean [0/Fe]-[Fe/H] trend of thick-disk stars with 
that of thin-disk members at [Fe/H] > —0.2. Nevertheless, we find that the kinematic membership 
criterion fails at separating perfectly the stars in the [O/Fe]- [Fe/H] plane, even when a very restric- 
tive kinematic separation is employed. Stars with "intermediate" kinematics (Pi < 0.7, P2 < 0.7) do 
not all populate the region of the [0/Fe]-[Fe/H] plane intermediate between the mean thin-disk and 
thick-disk trends, but their distribution is not necessarily bimodal. Halo stars (P3 > 0.5) show a large 
star-to-star scatter in [0/Fe]-[Fe/H], but most of it is due to stars with Galactocentric rotational 
velocity V < -200 kms" 1 ; halo stars with V > -200 kms" 1 follow an [0/Fe]-[Fe/H] relation with 
almost no star-to-star scatter. Early mergers with satellite galaxies explain most of our observations, 
but the significant fraction of disk stars with "ambiguous" kinematics and abundances suggests that 
scattering by molecular clouds and radial migration have both played an important role in determining 
the kinematic and chemical properties of solar neighborhood stars. 

Subject headings: stars: abundances — stars: atmospheres — stars: fundamental parameters - 
Galaxy: disk — Galaxy: evolution — Galaxy: formation 



1. INTRODUCTION 

Stars in the solar neighborhood exhibit a variety of 
kinematic and elemental abundance properties. Even 
though a majority of these stars, including the Sun, ap- 
pear to belong to a common stellar population, namely 
the Galactic thin disk, a significant number of nearby 
stars have been associated with the thick disk, whose 
members tend to be slightly more metal-poor and 
have enhanced a-element abundances (relative to their 
iron content). Both disk components are rotationally- 
supported, but the thick disk lags the thin disk, and 
it has a larger velocity dispersion. A small fraction of 
nearby stars having very old ages and low metallicities 
belong to the halo of the Milky Way. Even though re- 
stricted to a small region of the Galaxy, disentangling 
these local populations and studying in detail their prop- 
erties allow us to investigate the formation and evolution 
of the Milky Way galaxy. 

The obvious advantage of studying nearby stars is the 
fact that their analysis can be based on very high quality 
astrometric, photometric, and spectroscopic data. Un- 
fortunately, currently available observational resources 
allow us only to perform restricted surveys of nearby 
stars (distances closer than about 250 pc) at high spec- 
tral resolution (R = A/ A A > 20 000), implying severe 
sample biases and uncertain selection functions. The 
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latter prevent a straightforward comparison of observed 
chemical abundance patterns with models of Galactic 
evolution. This observational deficiency will be re- 
solved in the foreseeable future with ongoing surveys such 
as the Apache Point Obser vatory Galactic Evolution 
Experiment (AP OGE E, e.g.. lAllende Prieto et all 12051: 
iMaiewski et alj 120101: lEisenstein et al.l 120 111) , which is 
obtaining R ~ 20 000 infrared spectra of about 10 5 gi- 
ant stars, the Gaia- ESO Public Spectroscopic Survey 
(jGilmore et al.l [20121 ) . whose goal is to acquire visible 
R > 15 000 spectra of a similar number of carefully 
selected dwarf and giant stars in the Milky Way, and 
planned surveys such as the HERMES project (High 
Efficiency and Resolution Multi-Eleme nt Spectrograp h 
for the Anglo- Australian Telescope; e.g.. lFree"m an 2010), 
which will obtain R ~ 30 000 visible spectra of about 10 6 
stars. 

Analyses of nearly complete, unbiased data sets for 
nearby sta rs such as the Geneva-Copenhagen Survey 
CGCS. e.g-INordstrom et ~aT1[200l IH olmberg et all llM 
ICasagrande et al.ll2011h and low resolution spectroscopic 
studies of large samples (> 350 000) of faint and more 
distant stars such as the Sloan Extension for Galac- 
tic Understanding and Exploration (SEGUE/SPSS, e.g. 
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IpOOa Ide Jong et al.l [2010t iLee et ail W\% 
I2012al lbl: iSchlesinger et al.l l2012t l" are im- 
proving significantly our knowledge of Galactic struc- 
ture, dynamics, and chemistry, thus providing strong 
constraints to Galaxy formation and evolution scenarios. 
More details can in principle be seen in higher quality 
data, implying that it is still helpful at this moment to 
investigate chemical abundance patterns of few to several 
hundreds of nearby stars using high spectral resolution, 
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high signal-to-noise ratio data. 

Si nce pionee ri ng wo rk by, e.g.. iGilmore fc Re id (1983) 
and iSoubiranl (|1993l) . the picture of a Galactic disk 
composed of a kinematically cold thin disk and a kine- 
matically warm thick disk, as discrete p opulations, has 
been generally acc epted (but see, e.g., INorris h Rvanl 
119911 : iNorrisI 11999) rl Chemical abundance patterns of 
kinematically selected samples of sol ar neighborhood 
stars appear t o sup p ort this idea ( e.g.. iFuhrmannl 119981: 
Bensbv et all 120051: IReddv et all 120061: IRami'rez et al.l 
20071) . as well as the analysis of stars in low to 
medium spectral resolution surveys of much larger vol- 
ume s such as the Ra dial Velocity Experimen t (RAVE, 
e.g.JVeltz et al.ll200ll and SEGUE/SDSS (e.g. JLee et atl 
l201lF In all these works, as mentioned before, thick- 
disk stars are found to be generally older and statis- 
tically more metal-poor than thin-disk members (e.g. , 
Fuhrmannll 19981: IGratton et al.ll2000t IBensbv et al.1 120051: 
Reddv et al.l 120061: lAllende Prieto et al.l I2006D . The 
chemical abundances suggest an enhancement of a- 
elements relative to iron for thick-disk stars, although the 
amount of enhancement depends on the stellar metallic- 
ity (e.g.. IBensbv et al.|[2005l: IReddv et al.ll2006D . In fact, 
kinematically-selected thin- and thick-disk stars seem to 
have about the same a-element to iron abundance ratios 
at [Fe/H] > 00 

In the case of local, nearby stars, sample biases can 
be problematic when interpreting results obtained from 
the analysis of high quality spectroscopic data. The so- 
lar neighborhood is dominated by thin-disk stars. Thus, 
when thick-disk star samples are constructed, a strong 
kinematic criterion is often applied to avoid thin-disk 
"contamination," leading to thick-disk star samples that 
may be more representative of some extreme kinematic 
limit of the thick disk distribution rather than the av- 
erage thick disk. In most cases, an equivalently strong 
kinematic criterion is applied to construct a compari- 
son sample of thin-disk stars, thus removing anything 
that could have an intermediate behavior, either in the 
kinematics or in the chemical abundances. Studying in 
detail (i.e., with high resolution, high signal-to-noise ra- 
tio spectra) objects that have so far been arbitrarily re- 
moved from thin/thick disk chemical abundance trends 
could help us draw a more complete pictur e of Galac- 
tic ch emical evolution. In fact, the work by Fuhrmann 
(2008), who uses the magnesium abundance as indicator 
of a-element enhancement, has already hinted at the im- 
port ance of "transition" ob jects in the Galactic disk (see 
also iFuhrmamil Il998l l200l . 

Historically, theories for the formation and evolution 
of the Galactic disk have been divided into two generic 
groups, namely bottom-up and top-down models (e.g., 
lMajewskilfl993l: iFreeman fc Bland-Hawthorn! 12002ft . In 
the former, a thick disk is formed by secular processes 
internal to the Galaxy, while in the latter, mergers with 
satellite galaxies are usually called upon to explain the 

4 Disks are often referred to as the kinematically cold compo- 
nents of galaxies, in contrast to the hot halo and (classical) bulge 
components. In this paper, we refer to the thick disk as the warm 
sub-structure of the cold disk, to differentiate it kinematically from 
the colder thin disk. 

5 In this work we use the standard definitions: [X/Y] = 
log(N x /N Y ) - log(Nyr/N Y )e, and A x = log(Nx/Nu) + 12, where 
JVx is the number density of nuclei of the element X. 



existence of the thick disk. Until recently, the differ- 
ent merger scenarios seemed to be the only ones ca- 
pable of reproducing in detail most of the data ob- 
served in the solar neighborhood (e.g.. lAbadi et al.ll2003t 
iBrook et al.ll2005h . However, bottom- up models in which 
radial mixing is an important component of the process 
of Galactic disk formation have been shown to explain 
these data as successfully as the merger models, and in 
some cases even better (e.g.. iSchonrich fe; Binnevl f2009b: 
lLoebman et"aHl20Tl iBird et al.ll2012l 7 Nevertheless, re- 
cent results from the SEGUE/SDSS collaboration show 
that these models may have difficulties when dealing 
with observations of t he Galaxy on a larger-scale (e.g., 
iSchlesinger et al]|2012l ). Also, radial mixing alone cannot 
account for the existence of counter-rotating thick disks, 
which have been observed in some externa l galaxies (e.g., 
lYoachim fc Dalcantonll20 05l). Moreover. iMinchev et al.l 
(|2012t ) have used A^-bodv simulations to argue that the 
thickening of the Galactic disk may not be attributed to 
radial mixing. Despite all observational and theoretical 
efforts, we do not yet have a fully consistent picture for 
the formation and evolution of the Milky Way's disk. 

In this work, we measure oxygen abundances of a large 
number of nearby FGK stars, mainly dwarfs and sub- 
giants, to provide additional clues to the formation and 
evolution of the Galactic disk, and to a lesser extent, 
that of the Milky Way's halo. This work is an extension 
of a previous study in whi ch we presented oxy gen abun- 
dances of 523 nearby stars (jRamirez et al.ll2~007l hereafter 
R07). We have increased our sample size by more than 
300 objects and have improved the determination of stel- 
lar parameters and chemical abundances, as described in 
Sections [2] and [3] The oxygen abundance patterns de- 
rived from our data are presented in Section |H along 
with a discussion relevant to the chemical evolution of 
the Galaxy. Our conclusions and final remarks are given 
in Section [5j 

2. SAMPLE AND BASIC DATA 
2.1. Spectroscopic Data 

We have analyzed 897 spectra of 825 stars, in addition 
to 10 solar (day-sky and asteroid) spectra. Most of these 
spectra were taken by us, but we also used data from pub- 
lic archives. Multiple spectra for a number of stars are 
available from more than one source. Instead of combin- 
ing them, which would not be trivial given the differences 
in spectral resolution and the particular way in which the 
continuum normalization was made, in addition to differ- 
ences in wavelength coverage, we analyzed each spectrum 
independently and averaged the final results for the stel- 
lar parameters and elemental abundances. All spectra 
used have high resolution (R > 45, 000) and high signal- 
to-noise ratio (S/N > 100). This minimizes the impact 
of line blends and allows an accurate determination of 
local continua, which are important to reduce the obser- 
vational errors. 

Our spectra have been acquired using four instru- 
ment/telescope combinations, which are listed here, 
along with their corresponding references, in order of im- 
portance, as quantified by the number of spectra taken 
from each source: TS2/McD (R. G. Tull Coude spec- 
trograph, 2.7 m Telescope at McDonald Observatory; 
iTull et all fl995lh HRS/HET (High Resolution Spec- 
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Table 1 

Sources of Spectroscopic Data 



Source 


Facility /Instrument 


R 


S/N 


Spectra 


B03 


VLT/UVES 


80 000 


300-500 


55 


R03/R06 


McD/TS2 


60 000 


100-400 


369 


AP04 1 


McD/TS2 


60 000 


150-600 


106 




ESO/FEROS 


45 000 






R07 


HET/HRS 


120 000 


200-300 


52 




McD/TS2 


60 000 


150-250 


26 


AP08 


HET/HRS 


120 000 


100-200 


86 


R09 


McD/TS2 


60 000 


150-300 


69 


TW 


McD/TS2 


60 000 


100-300 


144 



Re ferences. — IBaenulo et af] ( 120031 . B03). | Reddy et al.| 

(2003, R03), IReddv et al.l 112009. R06), lATende Prieto et alj 
(2004, AP04).IRami'rez et al. (2007, R07), Allc nde PrietcT ^OOS. 
AP08). lRlmiirez et al.1 (j2009l. R09). This Work (TW). 

1 The number of spectra from AP04 correspond to the combined 
McD/TS2 and ESO/FEROS data set. 

trograph, 9.2 m Hobby-Eberly Telescope; iTullI I1998D . 
UVES/VLT (UV - Visual Echelle Spe ctrograph, 8 m Very 
Large Telescope: iDekker et alJEJoOO:. and FEROS/ESO 
(Fiber-feb Extend ed Range Optical Spect rograph, ESO 
1.52-m Telescope; iKaufer fe PasquinU 119981 ) . Details on 
the spectroscopic data reduction can be found in the ref- 
erences given below. 

The so urces of our spectra are li s ted in Table [TJ Spec- 
tra fro m IReddv et all (120031 12006D lAllende Prieto etaTI 
(|2004h . and IBaenulo et al.1 (|2003h were already ana- 
lyzed in R07, where we also included 78 stars from our 
own observations. For this work, we add the obs erva- 
tions of solar analog stars by lAllende Prietol (|2008l) and 
I Ramirez et al.1 (|2009ft . All these spectra have been re- 
analyzed in this work. In addition, we analyze new data 
we have obtained for 144 stars with the Tull spectrograph 
at McDonald Observatory. Most of these stars were se- 
lected to increase the number of stars with kinematic 
properties intermediate between those of the Galactic 
thin and thick disks. The reduction of these new data 
followed the exact same procedure described in R07. 

A few sample spectra are plotted in Figure [1] The 
strongest features seen in the solar spectrum are identi- 
fied. Besides showing the high quality of the data used, 
this plot roughly traces the behavior of the 777 nm Ol 
triplet lines, which we use as oxygen abundance indica- 
tors, across the region of stellar atmospheric parameter 
space covered by our sample. 

In all cases shown in Figure [TJ the 777 nm O I triplet 
is clearly identified and easily measured. At the spectral 
resolution of most of our data, line asymmetries due to 
surface convection are negligible, which allows us to use 
Gaussian or Voigt profile fits to accurately quantify line 
strengths. Line blending of the triplet is negligible in 
most FGK dwarf and subgiant stars, and it is important 
only for the coolest stars of our sample. Contamination 
by CN features can be seen in the spectrum of the cool 
giant star HIP 22489, as well as clear evidence for a blend 
in the middle line of the triplet, which is expected to have 
an intensity intermediate between those of the red and 
blue com ponents. This blend is probably due to a n Fe I 
line (e.g.. lTakeda et al.|[l998t iSchuler et al.ll2006aD . The 
middle line of the triplet was excluded from our cool K- 
dwarf and cool giant star analysis. 

It is important to note that there are no strong telluric 
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Figure 1. Sample spectra used in this work. Only a very small 
wavelength window centered at the O I 777 nm triplet is shown. 
The numbers in parenthesis are the atmospheric parameters T e g 
(in K), logg, and [Fe/H] of each star. 

features contaminating the wavelength region shown in 
Figure [TJ They would appear obvious as narrow spec- 
tral lines located at different wavelengths in each spec- 
trum, given that these sample spectra are corrected for 
the stellar radial velocities. Visual inspection of several 
telluric standard star observations confirms this state- 
ment. Moreover, the continuum normalization of this 
spectral window is simple since no extremely strong lines 
are present and relatively few spectral lines are visible at 
all, except in the case of the cool giant star HIP 22489. 

2.2. Photometry 

Visual magnitudes V and (B — V) colors were ex- 
tracted mainly fro m the General Catalogu e of Photomet- 
ric Data (GCPD, [Mermilliod et al1[l997l )FI For the few 
stars without Johnson's photom etry in the GCPD, we 
searched the Hipparcos catalog (jPerrvman et al.l I1997D 
for V and (B — V) values observed from the ground. 
These are values compiled from the literature by the 
Hipparcos team and not values inferred from transfor- 
mation equations using Tycho photometry, which are not 
so accurate. Cousins RI photometry was also extracted 
from the GCPD, where available. For stars fainter than 
V ~ 6, we used JHK S p hotometry from th e Two Micron 
All Sky Survey (2MASS. ICutri et al.ll200l . For brighter 
stars, 2MASS data are less reliable due to saturation. 

6 Available online at http://obswww.unige.ch/gcpd/gcpd.html 
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Stromgren (b — y) colors were also used, where available, 
as li sted in the GCPD and the Geneva-Copenhagen Sur- 
vey (jNordstrom et all 120041: ICasagrande et al.ll2011l) . 

2.3. Distance, Proper Motion, and Radial Velocity 

Most of our sample stars are included in the Hippar- 
cos catalog. We used the p arallaxes and proper motions 
from the new reduction by Ivan Leeuwenl (|2007l ). Stellar 
distances were derived from these trigonometric paral- 
laxes. For HD 144070, which is not included in the Hip- 
parcos catalo g, we used the ground-b ased measurement 
listed in the Ivan Altena et all |l995) compilation. The 
stars HIP 7751 and HIP 55288 are each known to be in 
wide visual binary systems. Their secondaries, for which 
we also have spectroscopic data and have therefore been 
analyzed, are fainter and not listed in the Hipparcos cat- 
alog. In these cases we adopted the parallaxes of the 
primaries for the secondaries. Our spectroscopic analy- 
sis confirmed the true binary nature of these two pairs; 
the radial velocities and elemental abundances we infer 
are in agreement for the two stars in each pair. 

Radial velocities were adopted from published catalogs 
and our own measurements. In R07, we provided radial 
velocities for a large fraction of our sample stars. For 
these objects, the radial velocities adopted here are from 
that paper. For the other stars, we used the average of 
values fr om the follow i ng sou rces: iBarbier-Brossat et al.l 
(|1994fl . iDuflot et ail (11995ft. t h e Geneva-Copenhage n 
Survey (INordstrom et al.l 120041). iRarmrez et al.1 (|2009h . 
IJenkins et al.1 (|2011[ ). lLatham et al.l (2002) , and our own 
measurements. Other sources for few stars, which 
were not found in a ny of the work s liste d ab ove, are 
iFamaev et all (|2005T) .IE vans fc Wild! (fl969l) . and lWilsonl 
(|1953D . 

Most stars have radial velocities with uncertainties un- 
der lkms^ 1 . A few are known spectroscopic binaries 
where the spectrum of the secondary is not detected, 
but important radial velocity variations are measured. 
For these objects we searched for orbital solutions in the 
lLatham et "all (|2002h survey and adopted their mean sys- 
tem velocity as the radial velocity of the star. A few 
of them do not have orbital solutions published and we 
therefore adopted a large error corresponding to the max- 
imum difference in radial velocity observed. 

2.4. Galactic Space Velocities and Kinematic 
Membership Criteria 

With the radial velocities and proper motions we de- 
rived heliocentric G alactic space velocities U, V, W fol- 
lowing the recipe by Uohnson fc Soderbloml (|1987l ). Er- 
rors in the input data were propagated as suggested 
by them. We note, however, that this recipe does 
not take into account the correlations between uncer- 
tainties in the astrometric quantities that are provided 
in the Hipparcos catalog. An impro ved error treat- 
ment is giv en in the Hipparcos c a talog ([Perrvman et all 
119971) or in lAllende Prieto et al.l (|2004fl . For stars that 
do not have an Hipparcos parallax we adopted the 
U, V, W values from the G eneva Copenhagen Survey 
(e.g.. ICasagrande et al.ll20"Tll ) and assumed conservative 
U, V, W errors of 2.5 km s _1 , as these values are not given 
for individual stars in the GCS. 

Membership probabilities (i.e., the probabilities of a 
star to be a thin disk, thick disk, and halo mem- 
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Figure 2. Distribution of visual magnitude and distance for our 
sample stars. Most of them are brighter than V = 10 mag and 
closer than d = 200 pc. 

ber) were computed as in R07 (their Section 3.3; see 
also iMishenina et al.1 [2004). In summary, stars of a 
given population are assumed to have a Gaussian dis- 
tribution in their Galactic space velocities, with mean 
U,V,W values a nd ve locity dispersions (a) adopted from 
iSoubiran et al.l (120031) for the thin and thick disks, and 
from lChiba & Beersl (|2000D for the haloQ The probabil- 
ity that a star with Galactic space velocity components 
U, V, W belongs to the thin disk (Pi), thick disk (P 2 ), or 
halo (P3) is thus given by: 

^ Ci 

P: 



(27r) 3 / 2 CTr/ i &Wi 
~U 2 



exp 



-0.5 



'V, 



W 2 



where 



(2) 



is a normalization constant which ensures that P^ = 1. 
The pi values are the relative number densities of thin- 
disk, thick-disk, and halo stars in the solar neighborhood, 
for which we adopt pi = 0.90, P2 = 0.08, and pz = 0.02. 

In this work, we associate stars with Pi > 0.5 with the 
thin disk (i.e., those with 50% or greater probability of 
being thin-disk members), P2 > 0.5 with the thick disk, 
and P3 > 0.5 with the halo. In order to compare our 
results with other works, in certain parts of this paper 
we will adopt, in addition to the former conditions, a 
"strong" kinematic criterion of P1/P2 > 1 for the thin 
disk and P2 /P1 > 10 for the thick disk (cf. iBensbv et al.l 
I2004L 120051 ). Note, however, that the latter criterion im- 
plies that a number of stars will be rejected from the 
analysis, which, as we have argued in Section [IJ is some- 
thing that could produce biased results. 

We provide in Table [5] visual magnitudes, trigonomet- 
ric parallaxes, radial velocities, Galactic space velocities, 
and membership probabilities along with their formal er- 
rors, for all our sample stars. Figure [2] shows the dis- 
tribution of apparent brightness and distance for these 
objects. A Toomre diagram is shown in Figure [H 

3. STELLAR PARAMETERS AND ABUNDANCES 

The majority of our sample objects are dwarf, turn- 
off, and subgiant stars with [Fe/H] > — 1.5 and a reliable 
trigonometric parallax measurement. The stellar param- 
eters for these stars, hereafter referred to as the "main" 

7 Thin disk: (Vi, aUi, &V\, aWi) = (-12, 39, 20, 20) kmr 1 , 
thick disk: (V 2 , <rU 2 , aV 2 , cW 2 ) = (-51, 63, 39, 39) km s" 1 , and 
halo: (V 3 ,<7U 3 ,crV3,crW3) = (-199, 141, 106, 94) kms" 1 . 
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Table 2 

Basic and Kinematic Data 1 



Star 


TJ 

•'mag 


7r (mas) 


V r (kms -1 ) 


U (kms -1 ) 


V (kms -1 ) 


W (kms -1 ) 


Pi 


P 2 


P 3 


HD 32071 


8.93 




8.3 ±0.5 


-24.0 ±2.5 


-25.0 ±2.5 


56.0 ±2.5 





77 





22 


0.00 


HD 59490 


8.68 




91.0 ±0.5 


-83.0 ±2.5 


-49.0 ±2.5 


-17.0 ±2.5 





70 





30 


0.00 


HD 67163 


8.05 




65.3 ±0.5 


-16.0 ±2.5 


-64.0 ±2.5 


-4.0 ±2.5 





70 





30 


0.00 


HD 82960 


8.53 




46.0 ±0.5 


28.0 ±2.5 


-58.0 ±2.5 


5.0 ±2.5 





80 





19 


0.00 


HD 130047 


8.56 




15.9 ±0.5 


41.0 ±2.5 


-63.0 ±2.5 


-21.0 ±2.5 





57 





43 


0.00 


HD 144070 


4.84 


43.0 ±4.2 


-32.6 ± 1.1 


-30.0 ±2.5 


-6.0 ±2.5 


-14.0 ±2.5 





99 





01 


0.00 


HD 170058 


8.42 




-21.2 ±0.5 


3.0 ±2.5 


-49.0 ±2.5 


-28.0 ±2.5 





86 





14 


0.00 


HD 171029 


8.24 




-37.8 ±0.5 


-22.0 ±2.5 


-65.0 ±2.5 


-7.0 ±2.5 





65 





34 


0.00 


HD 183490 


8.22 




-66.1 ±0.7 


-5.0 ±2.5 


-71.0 ±2.5 


0.0 ±2.5 





50 





19 


0.00 


HD 213746 


8.55 




-59.1 ± 1.9 


-2.0 ±2.5 


-71.0 ±2.5 


-15.0 ±2.5 





45 





54 


0.01 


HD 223723 


8.60 




4.3 ±0.5 


-85.0 ±2.5 


-35.0 ±2.5 


-26.0 ±2.5 





83 





17 


0.00 


HIP 171 


5.74 


82.17 ±2.23 


-36.9 ±0.5 


-8.0 ±0.7 


-73.1 ± 1.7 


-31.1 ± 2.0 





23 





76 


0.01 


HIP 348 


8.62 


17.43 ± 0.90 


18.1 ±0.5 


-94.0 ±6.1 


-12.1 ± 1.6 


-12.2 ±0.4 





91 





06 


0.00 


HIP 394 


6.11 


25.52 ±3.28 


4.6 ±0.5 


-124.8 ± 19.4 


-51.3 ±8.7 


-10.9 ± 1.5 





27 





71 


0.02 


HIP 475 


8.22 


19.01 ±0.82 


-28.4 ±0.5 


6.8 ±0.4 


-31.4 ±0.6 


-36.7 ± 2.0 





93 





07 


0.00 


HIP 493 


7.45 


26.93 ±0.56 


-45.6 ±0.5 


45.8 ± 1.1 


-34.2 ±0.4 


17.0 ±0.6 





95 





05 


0.00 


HIP 522 


5.70 


38.89 ± 0.37 


5.3 ±0.5 


-57.4 ± 1.0 


-37.5 ±0.6 


-15.3 ±0.5 





93 





07 


0.00 


HIP 530 


8.35 


9.38 ±0.71 


-31.8 ±0.5 


-28.2 ± 4.5 


-48.5 ± 2.4 


2.8 ±0.9 





92 





08 


0.00 


HIP 544 


6.10 


73.15 ±0.56 


-6.9 ±0.5 


-15.0 ±0.2 


-21.6 ±0.4 


-10.1 ±0.3 





99 





01 


0.00 


HIP 656 


8.15 


9.14 ±0.98 


-42.5 ±0.5 


-5.8 ±2.5 


-55.7 ± 1.9 


-29.0 ±3.6 





74 





25 


0.00 



1 Tabic [2] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 
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Figure 3. Toomre diagram for our sample stars. Error bars are 
shown only for stars with large velocity errors (> 20kms in the 
total Galactic space speed). Dotted lines show curves of constant 
speed, i.e., (U 2 + V 2 + W 2 ) 1 / 2 = const. 

sample, were determined as described in Sects. 13.11 to 
13.31 For the other stars, we used alternative methods, 
as explained in Section 13.41 Stellar parameter determi- 
nation is an iterative procedure because the end prod- 
ucts T e g, log g, [Fe/H] are inter-dependent. The results 
presented below were obtained after several iterations, 
which ensures that, for each star, the three parameters 
are self-consistent. An HR diagram of the full sample, 
using our final atmospheric parameters, is shown in Fig- 
ure m 

3.1. Effective Temperature 

Using as many as available of the (B — V), (V — Rc), 
(V - 7 C ), (Rc ~ Jc), (V - J), (V - H), (V - K s ), 
(J — K s ), and (b — y) colors, we derived the star's effec- 
tive temperatures f rom the metallicity-depe ndent color- 
ing calibrations by iCasagrande et al.l (|20100 . Since one 
T e fj value is obtained from each color, we computed a 
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Figure 4. HR diagram of our sample. The error bar on the top 
left shows the average size of our internal uncertainties, excluding 
systematic errors. 

weighted mean T c ff value and error for each star. The 
weights were computed using the photometric errors and 
th e mean accuracy o f eac h color calibration, as given 
bv lCas agrande ~et al.l (|2010T ) . Since the color-calibrations 
used depend on [Fe/H], the error in [Fe/H] that we obtain 
in Section [373] was also taken into account when estimat- 
ing the error in T c g. Most of our sample stars have at 
least five colors available, and their T e g values have an 
average 1 a color- to-color scatter of 48 K. 

Interstellar reddening is important only for stars more 
distant than about 100 pc, particularly if they are lo- 
cated close to the Ga lactic plane. As shown by, e.g., 
ILallement et al.l (|2003f ) . the Sun is located inside a "Lo- 
cal Bubble" of low interstellar gas density which ap- 
pears to be dust free. This region has a radius of about 
75 pc , but is far from sp herically symmetric (e.g.. iLerovl 
Il999t iWelsh et"aTl 120101 ) . We employed several inter- 
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stellar reddening maps to derive E(B — V) values for 
stars more distant than 75 p c (as described in detail in 
Ramirez & Melcndez 2005a, their Section 3.1). For stars 
closer than 75 pc we assumed E(B — V) = 0. 

The tempe r ature- color calibrations by 

ICasagrande et all (|2010| ) are based on effective tem- 
peratures obtained in the most recent implementation 
of the so-called infrared flux method (IRFM , e.g., 
iBlackwell fc Shallisl [19771 IBlackwell et all Il979t) . The 
zero point of the ICasagrande et all (|2010l ) T e s scale has 
been carefully calibrated and successfully tested using 
samples of well-known stars and precise spectrophoto- 
metric data. It surpasses previous implementat i ons o f 
the IRFM, in particular t hose b y lAlonso et ail (| 19961 ) 
and iRamfrez fc Melendezl (|2005aD . The latter was used 
in R07, and this constitutes one of the major improve- 
ments made to our work, particularly considering the 
fact that the zero point of the IRFM T e g scale was 
shifted by about 100 K in the more recent work by 
ICasagrande et al.l (|2010l ). This T c g offset translates into 
important changes to the [Fe/H] and [O/Fe] elemental 
abundance scales, as will be discussed later. 

3.2. Surface Gravity, Mass, and Age 

Surface gravities (logg), stellar masses (M), and ages 
(r) were computed by comparing the location of our 
sample stars in the HR diagram (T c g vs. My) with 
theoretical predictions based on stellar evolution calcu- 
lations. The absolute magnitudes were obtained from 
the stars' apparent magnitudes and trigonometric paral- 
laxes. We used the Yonsei-Yale isochrones, which take 
into account q-element enhancement a t low metallicities 
(|Yi et alll200Tl l2003HKim et alll2002D . 

In our age determination algorithm, each isochronc 
point is represented by a set of parameters (x,t,m, /), 
where x could be the surface gravity, mass, or age; t is the 
effective temperature, m the absolute magnitude, and / 
the iron abundance. To obtain the most likely parameter 
x given a set of observed quantities (T e ff, My, [Fe/H]) 
and their associated errors (AT c ff, AMy, A[Fe/H1), we 
use a probabilistic approach (see also, e.g.. iReddv et all 
20031 lAllende Prieto et all l200l iBaumann et al l 120101: 
Chaname fc Ramirez! 1 2012 : iMelendez et al.l 120121 . The 
probability that a given isochrone point corresponds to 
the observation is computed as: 

p(x, t, to, /) <x exp[-(T cff - tf / (V2AT cS ) 2 } x 

exp [-(My - m) 2 /(V2AM v ) 2 } x (3) 
exp[-([Fe/H]-/) 2 /(V2A[Fe/H]) 2 ] . 

Then, we calculate a probability distribution for the pa- 
rameter x as follows: 

x — x+8x 

P(x)= P(x',t,mJ), (4) 

x' —x — Sx 

where Sx is a small step for the parameter x while the 
sum extends also to all values of t, to, and / within a 
radius of 3 times the observational errors from the stel- 
lar parameters. Including isochrone points beyond these 
limits has no significant impact on the final result. The 
6x values correspond to a twentieth part (1/20) of the 
range of x available in each calculation. However, in 
cases where this Sx is too small considering the sampling 



of isochrones, we adopted 0.2 Gyr and 0.005 M for the 
age and mass distributions, respectively. The peak of the 
P(x) function, which corresponds to the x value of the 
bin in which P(x) is maximum, is adopted as the most 
probable surface gravity, mass, or age, while the 68 % 
and 96 % confidence limits are adopted, respectively, as 
the 1 a and 2 a Gaussian-like lower and upper limits. 

The isochrone grid used has very fine spacing, in 
particular in [Fe/H], where it is only 0.02 dex, which 
was achieved by employing the interpolation routine by 
IKim et all (|2002l ) . This allows us to determine with more 
precision the shapes of the age probability distributions 
without having to increase arbitrarily the error bar in 
[Fe/H], as it is sometimes done. The latter could intro- 
duce biases in the age determination scheme. 

A number of studies have shown that the probabilis- 
tic age determination approach that we have adopted is 
affected by b iases that can often overestimate the stel- 
lar ag es (e.g.. lPont fc Everll2004l : lJ0rgensen fc Lindegrenl 
2005). These biases are relevant for works where the ab- 
solute values of the ages are important, such as those 
dealing with the age-metallicity relation. Relative ages 
or sample chronology (i.e., simply sorting the stars ac- 
cording to their age) are less affected by these biases, 
yet not completely free from them. In any case, we 
note that a comparison of our ages with those deter- 
mined using Bayesian approaches that at tempt to min- 
imize those biases, specifically those by Ida Silva et all 
(2006) and Casagrand e et all (1201 ID. reveals surpris- 
ingly g ood agreement, as shown bv lChaname fc Ramirez! 
(2012, their Figure 7), although in part this could also 
be due to the use of different isochrone sets which may 
compensate the differences due to the statistical biases. 

3.3. Iron Abundance 

We determined [Fe/H] values using a standard spec- 
troscopic approach. The 2010 version of the spectrum 
synthesis code MO OG (ISnedenl \WfW\ an d MARCS 
model atmospheres (|Gustafsson et all 120101 ) with stan- 
dard composition (i.e., including a-element enhancement 
at [Fe/H] < 0)0 interpolated linearly to the atmospheric 
parameters of each star, were used. A large number of 
neutral (Fei, up to 128) and singly-ionized (Fen, up to 
16) iron lines were used for the determination of [Fe/H]. 
The version of the spectrum synthesis code and the 
model atmosphere grid adopted are different from those 
used in our previous work. In R07, we used the 2002 
version of MOOG and Kurucz models with solar-scaled 
composition, i.e., without a-elemcnt enhancement. 

Equivalent widths were measured using an automated 
routine capable of identifying blends03 Spectral lines are 
assumed to have Gaussian shape, which is a poor approx- 
imation for strong lines with extended wings. However, 

8 http://www.as.utexas.edu/~chris/moog.html 

9 Downloaded from http://marcs.astro.uu.se These standard 
co mposition MARCS m odels use scaled solar abundances (those 
by IGrevesse et al. 2007), but with [a/Fe] and [O/Fe] abundance 
ratios increasing linearly from +0.00 to +0.40 as [Fe/H] decreases 
from +0.00 to -1.00. For the [Fe/H] < -1.00 models, [o/Fe] and 
[O/Fe] are set equal to +0.40. 

10 This routine, "getew_xy.pro," written in IDL, is available 
online at http://hebe.as.utexas.edu/stools Also available at this 
website is the interpolation routine tor MARCS model atmospheres 
used in this work, "mmod.pro," which was adapted from the code 
"kmod.pro" that works with Kurucz models. 
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we use only relatively weak lines for our iron abundance 
analysis, thus minimizing this systematic error. Only for 
very metal-poor stars, the equivalent widths were mea- 
sured manually using IRAF's task splot Fl because even 
small errors in the determination of the local continuum 
result in a large percent error in the measured equivalent 
width of very weak lines. For very metal-poor stars we 
find that the local continuum of weak lines is best deter- 
mined visually (rather than assuming that it is exactly 
at 1, which would imply a perfect continuum normaliza- 
tion). Since fewer lines are available for these objects, the 
higher precision in the equivalent width measurements 
that can be easily achieved with this approach results in 
more reliable average elemental abundances. 

Average [Fe/H] values were obtained from all Fel and 
Fe II lines analyzed, not as the average of Fe I and Fe n 
abundances measured separately, but giving equal weight 
to each of the iron lines (both neutral and singly ionized) 
measured. This procedure is safe because the mean iron 
abundances derived from Fel and Fen separately agree 
reasonably well, as discussed in Section [3.3.3l The [Fe/H] 
values used for tracing iron versus oxygen abundance 
trends in most of Section HJ however, are those inferred 
from the Fen lines, because they are less affected by sys- 
tematic errors on both the spectral line calculations and 
on the adop ted stellar parameters, as will be discussed in 
Section [4. II Errors in [Fe/H] were determined by adding 
in quadrature the standard error of the line-by-line iron 
abundance scatter and the small changes in [Fe/H] due 
to ±1 a variations of T e g and log g. 

3.3.1. Line Selection and Solar Analysis 

Our iron line selection and atomic data adopted are 
largely from R07. A few additional spectral lines 
and some minor updates to the atomic data were 
made following the works by iMelendez et all (|2009() and 
lAsplund et aU ([20091 ). The former is a high precision 
elemental abundance study of solar twin stars whereas 
the latter represents the most complete solar abundance 
analysis to date. The addition of these lines allowed us 
to have a better coverage of excitation potential and line 
strength with very clean features. The pr essure broaden- 
ing da mpin g constants adopted are from B arklc m~et al.l 
( 2000) and lBarklem fc Aspelund- Johansson! (|2005ft . but 
for lines without damping constants computed by them, 
we used the values obtained from Unsold's formula, mul- 
tiplied by a factor of 3. Our adopted iron line list is given 
in Tabled 

Ten spectra were available for a solar analysis. Half 
of them correspond to daysky observations whereas the 
other half are reflected Sun-light observations of bright 
asteroids (Ceres, Iris, Pallas, and Vesta). All these spec- 
tra are of very high quality and the measured equivalent 
widths agree very well, with the exception of one daysky 
spectrum, as discussed below. The average standard de- 
viation of EW values measured for the same line in all 
solar spectra is 1.8 mA, which corresponds to a standard 
error for the mean EW of 0.6 mA. Since we used our au- 
tomated routine to measure the solar EWs, this result 

11 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation - http://iraf.noao.edu 



Table 3 



Iron Line List 1 


A (A) 


Species 


EP (eV) 


logs/ 


EW Q 


4630.12 


Fel 


2.280 


-2.52 


72.4 


4635.85 


Fel 


2.850 


-2.34 


54.1 


4683.56 


Fel 


2.830 


-2.41 


55.9 


4690.14 


Fel 


3.690 


-1.61 


55.5 


4745.80 


Fel 


3.650 


-1.27 


78.3 


4576.33 


Fen 


2.844 


-2.95 


64.4 


4620.51 


Fen 


2.828 


-3.21 


52.4 


4629.34 


Fen 


2.810 


-2.28 


97.2 


5197.57 


Fen 


3.231 


-2.22 


81.7 


5234.62 


Fen 


3.221 


-2.18 


83.9 



1 Table[3]is published in its entirety in the elec- 
tronic edition of the Astrophysical Journal. A 
portion is shown here for guidance regarding its 
form and content. 

gives us a rough estimate of our automated EW mea- 
surement errors (about 1%). 

Only one of our solar spectra revealed a small but non- 
negligible offset in the measured EWs. This daysky spec- 
trum was taken with the HRS /HET at sunset and point- 
ing the telescope to the East, resulting in a large separa- 
tion angle between the area of the sky observed and the 
solar positio n. As shown by IGrav et al.l (|2000l . see also 
Section 2 in lRamfrez et al. l 1201 If ) . scattered light in the 
Earth's atmosphere can affect significantly the shapes 
and also the strengths of spectral lines observed under 
these conditions (i.e., at large angular separations be- 
tween the Sun and the area of the sky observed) Q Thus, 
before averaging the solar EWs, we corrected those of 
this spectrum by increasing them by 3 %, which brought 
them to close agreement (on the average) with the other 
solar spectra. We note also that this correction factor 
is in reasonably g ood agreement with that suggested by 
IGrav et ail (|2000j) . After applying this correction, for 
each spectral line we averaged the EWs measured in all 
10 spectra and adopted the mean value as the solar ref- 
erence EWs, which are listed in column 5 of Table [3] 

Using our iron line list and a MARCS model atmo- 
sphere, we derive a solar iron abundance A-p c = 7.46 ± 
0.06 using only the Fe I lines, and Ap c = 7.44±0.08 using 
only Fe n lines. Our inferred solar iron abundance is thus 
in reasonably good agreement with estim ates found else- 
where , in particular with that derived by lAsplund et al.l 
(2009) using a three-dimensional (3D) time-dependent 
hydrodynamic solar photosphere simulation. Note that 
3D corrections to the ID analysis (which we employ) are 
very small for iron lines in the solar case. 

3.3.2. Mtcroturbulence 

As is common practice, the microturbulent veloc- 
ity (v t ) was obtained by removing the correlation be- 

12 Interestingly, the other daysky spectra have EWs that are 
in good agreement with those measured in the asteroid spectra, 
which are point sources observed during night time. These daysky 
observations were made with the Tull spectrograph by letting scat- 
tered sunlight pass through a "solar port" which points towards the 
zenith, and a few hours before sunset, implying zenith-Sun angles 
between about 40 and 70 degrees. Although the solar spectral lines 
are expected to be distorted under these conditions, the impact on 
the EW values is likely smaller. Note also that this effect is better 
seen at higher spectral resolution. The HET-HRS daysky spectrum 
has R = 120 000 whereas the Tull daysky spectra have R = 60 000. 
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Figure 5. Top panel: microturbulent velocity inferred from our 
analysis of Fel lines as a function of T e g. Middle panel: as in 
the top panel for logg. Bottom panel: comparison of microturbu- 
lent velocities, inferred from the Fel lines, and obtained from the 
linear regression (Equation [5J- The dotted line shows the 1-to-l 
correspondence. The stars excluded from the regression leading to 
Equation [5] are not shown in this figure. 



tween iron abundance and reduced equivalent width 
(log(EW/X)) of Fel lines. For the solar reference we 
used the absolute abundances Ap \ and obtained Vt = 
1.03kms _1 . For the other stars we used the relative iron 
abundances, i.e., the [Fe/H] values from the Fel line anal- 
ysis. This procedure is more reliable than using absolute 
abundances because errors in the atomic data, and to a 
lesser extent in the model atmospheres, are minimized 
with our line-by-line differential analysis. 

The resulting vt values correlate strongly with T c g and 
log 5, as shown in Figure [5] for the dwarf and subgiant 
stars. Giant stars (whose analysis is discussed later in 
Section I3.4.2[) are excluded from Figure [5] because they 
are at the cool T e g end but have large v t values (from 
about 0.8 to 1.6 kms -1 ) compared to the cool K-dwarfs 
(v t < 0.5 km s _1 ). The Vt values show only a weak depen- 
dence on [Fe/H]. The dependency of Vt on stellar param- 
eters that we obtain is consistent with the fact that the 
microturbulent velocit y scales wel l with the strength of 
surfa ce convection (cf. iGravl 119781 : iNordlund &; Dravinsl 
1990). We performed a linear regression for vt on these 
three parameters, using an iterative 2.5 a clipping algo- 
rithm to remove outliers. We found the following mean 



relation: 

v t (km s" 1 ) = 1.163 + 7.808 x 10" 4 (T cff - 5800 K) 

-0.494(log. 9 - 4.30) - 0.050 [Fe/H] , (5) 

with a 1 a error of 0.12kms _1 . The limits of applicability 
of this relation are: T cff = 4940-6750 K, logg = 3.6-4.7, 
and [Fe/H] = -1.2 to [Fe/H] = +0.4. 

A comparison of Vt values derived using Equation [5] 
with those obtained in the iron line analysis is shown in 
the bottom panel of Figure [5j For the stars that were 
removed from this regression due to the 2.5 a clipping 
constrain (about 15% of the sample), we adopted the Vt 
values inferred from Equation [5] rather than that derived 
from their Fe I line analysis. Note that v t is a parameter 
that could be severely affected by the few spectral lines 
on the extremes of the reduced equivalent width distri- 
bution and therefore the mean relation is more reliable 
for the outliers. The latter is also true for cool K-dwarfs 
(T c ff < 5000 K), for which their vt values were inferred 
mainly by extrapolating EquationJSJ Note, however, that 
in none of these cases we obtained an unphysical (i.e., 
negative) Vt- 

Compared to other vt parameterizations found in the 
literature, ours is in goo d qualitative agre e ment . For 
example, the for mulas bv lEdvardsson et al.l (|1993l their 
Equation 9) and lAllende Prieto et al.l (|2004L their Equa- 
tion 2) both have a positive slope for the T e g dependence, 
but a negative one for the logg de pendence. Quan- 
titativ ely, the T e ff dependence of the lEdvardsson et al.l 
(1993) formula for dwarf stars (log g = 4.5) is almost 
identical to ours (mean difference of only 0.03kms _1 
with a 1 a standard deviation of 0.01 kms -1 , with our 
Vt values being smaller). Th e logq dependenc e , how - 
ever, is stronger according to lEdvardsson et ail (|1993|) . 
who predict w t ~ 1.9 km s -1 at \ogg = 3.8 while our pa- 
rameterization sugge s ts Vt — 1.3kms -1 . The formula by 
lAllende Prieto eta l. (2004), on the other hand, implies 
even shallower T e ff and logg dependencies, but a log<? 
dependence that i s mor e similar to ours than that by 
lEdvardsson et al.l (|1993l ). 

3.3.3. Spectroscopic Equilibrium 

Ideally, the average iron abundance inferred from Fe I 
and Fell lines separately should agree. In addition, the 
Fe I abundances should not correlate with the excitation 
potential of the lines. In many studies, these conditions 
are used to fine-tune T e ff, logg, and Vt, given their sensi- 
tivity to those trends. We prefer to avoid that approach 
because it is strongly sensitive to model uncertainties 
whereas our atmospheric parameter determination uses 
techniques that are known to be less affected by those 
systematic errors. Fine-tuning by ionization/excitation 
balance, however, can be powerful for the analysis of stars 
which are very similar to the Sun, when the latter is used 
as reference star, because systematic errors are essen- 
tially cancelled out in a line-by-line differential analys is 
(|Melendez et al.ll2?M 12013 IRamfrez eTalMM . I2IT1I . 

The Fe 11 minus Fe 1 abundance difference from our iron 
abundance measurements is shown in Figure [6] This plot 
shows the stars from our "main" sample, but it does not 
include those for which vt was inferred from Eq. [5] (i.e., 
stars which were sigma-clipped from the fit are not plot- 
ted, which includes many with T e ff < 5000 K). There is 
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Figure 6. Test of iron ionization equilibrium. The difference 
between the mean iron abundance determined from Fell and Fcl 
lines is shown as a function of stellar parameters. 

significant scatter in Figure [6] but no obvious trends or 
offsets are observed. The average Feii-Fei difference is 
0.01 ± 0.06 dex. This is contrary to the case of R07, who 
found an overall offset of +0.06 dex for the Fen abun- 
dances relative to Fe I, as well as severe discrepancies for 
cool metal-rich dwarfs. 

A more detailed inspection of the trends shown in 
Figure [6j particularly adding the stars with v t derived 
from Eq.[5l which includes several cool K-dwarfs, reveals 
that there is a trend of higher Fe II abundances for the 
cool metal-rich dwarfs in our data, as in R07. In fact, 
the mean Fe n minus Fe I difference for dwarf stars with 
Tcff < 4800 K is +0.09 ±0.07, a number that increases to 
+0.15+0.05 if the sample selection is further restricted to 
stars with super-solar metallicity. This severe iron ion- 
ization problem is seen only in cool metal-rich K-dwarfs, 
and it has also been detected in other works, as will be 
discussed in detail in Section B~T1 

Excluding the stars mentioned above, we do not detect 
an overall offset in the iron abundances inferred from Fe I 
and Fe II lines separately at all T e g values, as we found 
in our previous work. This is most lik ely due to the use 
of the improved IRFM T c g scale by iCasagrande et al.l 
( 2010), who shifted the ze ro point by about +100 K wit h 
respect to the T e e scale bv lRamfrez fe Melended ((2005b). 
A +100 K offset in the T e g values used is almost the 
exact amount necessary to remove the Fe n minus Fe I 
abundance offset in the R07 data. 

Figure [7] shows the [Fe/H] vs. EP slope for Fei lines. 
Our Fe I linelist consists of features with EP from about 
to 5 eV. The EP slope can be severely affected by one or 
two unreliable lines on the extremes of the EP distribu- 
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Figure 7. Test of Fei excitation equilibrium. The slope of abun- 
dance versus excitation potential for Fe I lines is shown as a function 
of stellar parameters. 

tion, which most likely explains the outliers in Figure [71 
The bulk of the data, however, reveals clear trends with 
all three parameters. Although excitation balance (i.e., 
no EP slope) appears to be satisfied around the solar 
values, the EP slope is clearly positive for relative low 
log<? < 4.3 and negative for relative low [Fe/H] < —1.0. 
It is also low at cool temperatures (T e g < 5500 K). This 
implies that the T e g values inferred using the condition 
of excitation balance are warmer than the effective tem- 
peratures obtained with the IRFM color calibrations for 
l°gff ^ 4.3 stars, but cooler for [Fe/H] < —1.0 and/or 
T eff < 5500 K stars. 

Our atmospheric parameters do not satisfy perfectly 
the so-called spectroscopic equilibrium conditions, yet 
there has been a significant improvement with respect to 
R07 in that ionization balance is no longer an important 
issue, except at the cool metal-rich end. The discrepancy 
in the excitation equilibrium is likely due to systematic 
errors related to model atmospheres and departures from 
local thermodynamic equilibrium (LTE) in the spectral 
line formation calculations. Finding a solution to this 
problem is beyond the scope of this paper. We refer the 
rea der to the recent works bv iBcrgcma nn et al.l (2012) 
and lLind et al.l (J2012I ) for a detailed investigation of these 
issues. 

3.4. Special Cases 
3.4.1. No Trigonometric Surface Gravity 

A number of our sample stars (11) are not included 
in the Hipparcos catalog and do not have ground-based 
measurements of their trigonometric parallax. They can 
be found in the Tycho catalog, but their parallaxes have 
large errors ranging from 30 to 100 %, with most of them 
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showing an error in the parallax of at least 50%. Thus, 
we cannot use the technique described in Section [3T21 to 
derive their logg values with high precision. For this 
sub-sample of stars (the "noplx" sample) we obtained 
logg by forcing the mean Fei and Fen abundances to 
agree (ionization equilibrium). The effective tempera- 
tures were still obtained as described in Section [XT] 

We used the same approach for two stars (HIP 5134 
and HIP 61272) which, even though they have an entry 
in the Hipparcos catalog, inspection of ionization equi- 
librium (Fe I minus Fe II abundance difference) and the 
wings of strong lines such as Mgi6 suggest a severely 
erroneous parallax. 

3.4.2. Giants 

The lCasagrande et al.l ()2010[ ) color- T e s relations apply 
only to dwarf and subgiant stars. No extension of this 
T e g scale for giants has been published yet. Therefore, for 
a number of stars on the red giant branch (27, the "giant" 
sample) we determined T s using excitation equilibrium. 
However, we did not force the EP slope to be exactly 
zero but O.OldexeV -1 because this value is the average 
observed in the relatively low logg stars from the main 
sample (see Figure [7]) . In this way, the T e s values for 
giant stars are forced to be roughly on the same scale as 
that of the dwarf and subgiant stars. However, we cannot 
guarantee that this will be the case for all the giant stars, 
given that they cover a very wide range of radii, implying 
significantly different atmospheric structures compared 
to stars at the base of the red giant branch. Neverthe- 
less, we note that for Arcturus (HIP 69673), which is one 
of the stars with the lowest log g in our sample, we derive 
T e g = 4215 ± 90 K, a value that is in good agreement, 
within 1 er, with that obtai ned in the very detaile d anal- 
yses of this giant star bv |Griffin fe Lvnas-Gravl (1999. 
T off = 4290±30K) and lRarmrez fe Allende Prietol |201lL 
T off = 4286 ±30 K). 

3.4.3. Other 

A few stars (18) presented several different difficulties 
when attempting to measure their atmospheric parame- 
ters with the techniques described above. In some cases 
very few lines were available, making the determination 
of v t extremely difficult. These could be not only very 
metal-poor stars, but also objects for which our spectra 
are of relatively low quality. In other cases, the published 
parallaxes returned unphysical logg values, as confirmed 
from inspection of the gravity-sensitive Mg I b features at 
5180 A or from their odd location in the HR diagram. 
For these stars (the "manual" sample), whenever pos- 
sible, we constrained T e g and/or logg from excitation 
and ionization balance; otherwise we adopted the aver- 
ages of values found in the literature, as compiled by 
iSoubiran et "all (|2010l ). 

3.5. Catalog of Atmospheric Parameters 

Our final atmospheric parameters are listed in Table [4j 
Results from multiple spectra of the same star have al- 
ready been averaged in this table. The latter was applied 
to 68 objects, with an average spectrum-to-spectrum 1 a 
scatter of less than 0.02 dex in [Fe/H]. For logg and T e ff, 
this scatter is even smaller: 0.01 dex and 2K, respec- 
tively. This is due to the fact that logg and T e g are de- 
rived almost independently from the spectrum; [Fe/H] is 



an input quantity that has only a minor impact on their 
derivation. Photometric and astrometric errors dominate 
the uncertainties in T c g and logg, respectively. Since 
the spectrum-to-spectrum scatter of the derived stellar 
parameters is small compared to the internal error (see 
below), we conclude that the heterogeneity of our spec- 
troscopic data set has a negligible impact on our results. 

The median errors of the stellar parameters (internal 
only, i.e., not including estimates of systematic errors) 
are: AT cff = 44 K, A logg = 0.04, and A [Fe/H] = 0.05 
(the way in which the individual errors were determined 
are described in previous sections). 

In Table 0] a flag for vt tells whether the microtur- 
bulence was derived from the Fei analysis or from the 
linear regression (Equation [5]) . Similarly, we include a 
column with information on how exactly were the at- 
mospheric parameters derived (column "sample"). Our 
adopted st ellar masses and ages, derived as described 
in Section 13.21 are given in Table [5] For stars with 
no trigonometric parallax available, hence with unknown 
absolute magnitude, the exact same mathematical proce- 
dure explained in Section l3~2l was used, replacing My, the 
observed absolute magnitude, with the spectroscopically 
measured logg, and m, the absolute magnitude from the 
models, with the logg values from the isochrones. 

3.6. Oxygen Abundance 

As shown before, the Oi triplet at 777 nm is a 
very strong feature in solar-type stars, easily detectable 
in high resolution, high signal-to-noise spectra. The 
equivalent widths of these three lines can be mea- 
sured with high precision for most nearby FGK-type 
stars (cf. Figure [TJ). Accurate transition probabilities, 
based on quantum m e chani ca l calculations, are availabl e 
fro m iHibbert et al.1 (1199 lh. iButler fe Zeippenl (|1991ft . 
and lBiemont fe Zeippenl (|1992h . These published values 
are reasonably consistent with each other and we there- 
fore adopted their averages: logg/ = 0.352,0.223,0.002 
for A = 7771.9, 7774.2, 7775.4A, respectively. 

However, it has long been known that the oxy- 
gen infrared triplet lines form in conditions that are 
not well reproduced under the PTE assumption (e.g., 
lEriksson fe Toftlll979l: rKisehnanlll993l) . Statistical equi- 
librium calculations show that in the solar case, the 
LTE oxygen abundance inferred from the O I 777 nm 
triplet is too high by ~ 0.2 dex (~ 50%) relative to 
other oxy gen features less affected by model uncertain- 
ties ( e.g.. iKiselmanl 119931 : iTakedal Il994t iRarmrez et al.l 
120071: iFabbian et al.ll2009l ). Non-LTE effects are stronger 
for warmer and lower surface gravity stars, owing to 
the stronger radiation fields and lower densities, which 
make radiative transitions more important relative to 
collisional excitation. Nevertheless, even for the coolest 
dwarf stars in our sample, non-LTE effects can introduce 
(predicted) errors of almost 0.1 dex. On the other hand, 
the LTE oxygen abundances of the warmest stars in our 
sample are overestimated by about 0.5 dex, i.e., by more 
than a factor of 2. 

In R07, we used an oxygen model atom composed of 
54 levels pl us the continuum and 242 transitions, con- 
structed by I Allende Prieto et"aL1 (|2003l ) , along with the 
computer codes TLUS TY and SYNSPEC (|Hubenvil 19881: 
iHubenv fe Lanzlll995f) to calculate level populations for 
the oxygen atom and emergent flux synthetic O I 777 nm 
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Table 4 

Atmospheric Parameters 1 



Star 


T eff (K) 


lo£ [cetsI 


[Fe/Hl v f 

M- *-'/ - L - 1 -J v L 


(km s 


M cal. 2 [Fc/Hl„ T 


[Fe/Hl „ „ 

L / JFc II 


sample 


HD 32071 


5894 ± 38 


4.23 ±0.12 


-0.34 ±0.04 


1.35 


-0.34 ±0.04 


-0.33 ±0.07 


noplx 


HD 59490 


5627 ± 27 


4.49 ±0.11 


-0.17 ±0.04 


0.93 


-0.17 ±0.03 


-0.17 ±0.08 


noplx 


HD 67163 


5831 ± 29 


4.26 ±0.09 


0.03 ±0.03 


1.26 


0.03 ± 0.03 


0.03 ± 0.04 


noplx 


HD 82960 


6450 ± 129 


3.95 ±0.20 


0.04 ± 0.04 


1.85 


0.04 ± 0.04 


0.02 ± 0.07 


manual 


HD 130047 


5590 ± 33 


4.31 ±0.12 


-0.02 ± 0.04 


0.91 


-0.02 ± 0.04 


-0.02 ±0.09 


noplx 


HD 144070 


6532 ± 87 


4.16 ±0.08 


0.02 ±0.10 


1.81 


+ 0.01 ±0.09 


0.07 ±0.13 


main 


HD 170058 


5896 ± 46 


4.21 ±0.14 


0.13 ±0.05 


1.29 


0.13 ±0.04 


0.13 ± 0.10 


noplx 


HD 171029 


5562 ± 54 


3.85 ±0.21 


-0.36 ±0.07 


0.97 


-0.36 ±0.07 


-0.36 ±0.08 


noplx 


HD 183490 


5759 ± 43 


3.81 ±0.14 


0.06 ±0.05 


1.25 


0.06 ±0.05 


0.07 ± 0.05 


noplx 


tiD 213746 


5813 ± 49 


A i o 1 n 1 o 

4.12 ± 0.18 


0.05 ± 0.06 


1.34 


A fir 1 ri fir 

0.05 ± 0.05 


0.05 ± 0.11 


noplx 


HD 223723 


5946 ± 41 


4.28 ±0.10 


0.01 ±0.03 


1.28 


0.01 ±0.04 


-0.01 ±0.02 


noplx 


HIP 171 


5510 ±66 


4.46 ± 0.01 


-0.76 ±0.06 


0.96 


-0.76 ±0.05 


-0.80 ±0.11 


main 


HIP 348 


5746 ± 55 


4.38 ± 0.06 


-0.19 ±0.04 


1.13 


-0.20 ±0.03 


-0.18 ±0.07 


main 


HIP 394 


5635 ± 37 


3.78 ± 0.07 


-0.48 ± 0.04 


1.24 


-0.48 ± 0.03 


-0.47 ±0.06 


main 


HIP 475 


5836 ± 72 


4.35 ± 0.05 


-0.34 ±0.06 


1.18 


-0.33 ±0.05 


-0.41 ± 0.10 


main 


HIP 493 


5960 ± 44 


4.41 ± 0.03 


-0.20 ±0.04 


1.16 


-0.20 ±0.03 


-0.20 ± 0.07 


main 


HIP 522 


6251 ± 44 


4.21 ±0.02 


0.05 ± 0.05 


1.56 


+ 0.04 ± 0.05 


0.07 ± 0.06 


main 


HIP 530 


5866 ± 40 


3.90 ±0.05 


-0.01 ±0.05 


1.36 


-0.02 ± 0.03 


0.03 ± 0.09 


main 


HIP 544 


5458 ± 40 


4.52 ±0.02 


0.14 ±0.06 


0.78 


+ 0.13 ±0.04 


0.17 ±0.13 


main 


HIP 656 


5805 ± 39 


3.82 ±0.06 


-0.24 ±0.06 


1.32 


-0.25 ±0.05 


-0.15 ± 0.07 


main 


Table g] is 


published in its entirety in the 


electronic edition 


of the 


Astrophysical Journal. A portion is shown 


here for 



guidance regarding its form and content. 

A + sign is written in this column if the adopted microturbulcnt velocity for the star was obtained from the linear regression 
(Equation [5j. 



Table 5 

Age and Mass 1 



Star 


Age (Gyr) 


-la 


±1<T 


-2<r 


+2cr 


Mass (M Q ) 




-la 


+ 1<T 




-2a 


H 


-2a 


HD 32071 


9.09 


7.12 


9.68 


3.92 


10.75 


0. 


.942 


0. 


929 





.992 


0. 


.906 


1. 


.054 


HD 59490 


6.43 


3.44 


10.64 


1.15 


12.89 





.918 


0. 


895 





.936 





.877 





.956 


HD 67163 


7.31 


5.64 


7.74 


3.31 


8.46 


1. 


.032 


1. 


016 


1 


.062 


1. 


.000 


1 


.096 


HD 82960 


2.01 


1.71 


2.89 


1.06 


3.84 


1. 


.313 


1. 


313 


1 


.692 


1. 


.215 


1 


.873 


HD 130047 


11.65 


7.10 


12.17 


2.88 


13.39 





.936 


0. 


914 





.957 


0. 


.898 





.993 


HD 144070 


2.14 


1.73 


2.55 


1.12 


3.00 


1. 


.384 


1. 


334 


1 


.453 


1. 


.283 


1. 


.520 


HD 170058 


5.22 


4.01 


6.29 


1.81 


7.20 


1. 


.102 


1. 


078 


1 


.195 


1. 


.037 


1. 


.281 


HD 171029 


6.89 


4.37 


12.53 


2.24 


14.36 





.940 


0. 


906 


1 


.242 


0. 


.851 


1 


.535 


HD 183490 


3.90 


3.24 


6.60 


2.21 


8.36 


1. 


.242 


1. 


139 


1 


.446 


1. 


.052 


1. 


.637 


HD 213746 


7.16 


5.04 


8.03 


2.65 


9.20 


1 


.058 


1. 


013 


1 


.158 





.982 


1. 


.323 


HD 223723 


5.70 


3.96 


6.22 


1.83 


7.04 


1 


.078 


1. 


053 


1 


.119 


1. 


.031 


1. 


.169 


HIP 171 


14.46 


12.65 


14.56 


10.92 


14.78 





.759 


0. 


.751 





.771 





.744 





.787 


HIP 348 


8.67 


5.28 


10.20 


2.63 


12.19 





.930 


0. 


913 





.963 


0. 


.891 


0. 


.993 


HIP 394 


5.00 


3.93 


7.22 


3.01 


9.34 


1 


.132 


1. 


053 


1 


.257 





.977 


1. 


.369 


HIP 475 


8.75 


7.09 


11.47 


4.45 


13.32 


0. 


.923 


0. 


891 





.954 


0. 


.863 


0. 


.989 


HIP 493 


5.39 


3.21 


6.06 


1.78 


7.21 


1 


.004 


0. 


986 


1 


.034 


0. 


.966 


1. 


.057 


HIP 522 


3.00 


2.70 


3.40 


2.37 


3.63 


1. 


.244 


1. 


232 


1 


.258 


1. 


.222 


1. 


.274 


HIP 530 


4.90 


4.40 


5.60 


3.85 


6.16 


1 


.222 


1. 


184 


1 


.284 


1. 


.142 


1. 


.333 


HIP 544 


3.51 


1.59 


6.23 


0.46 


8.80 


0. 


.950 


0. 


919 





.976 


0. 


.892 


0. 


.995 


HIP 656 


4.19 


3.51 


5.39 


2.90 


6.51 


1 


.254 


1. 


178 


1 


.346 


1. 


.109 


1. 


.429 



1 Tablc[5]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. 



triplet lines in a grid of stellar atmospheres. Non-LTE 
corrections to the abundances derived using the LTE ap- 
proximation were then computed for this grid using a 
wide range of oxygen abundances. We used this grid of 
non-LTE corrections to interpolate linearly to the stellar 
parameters of our sample. The grid has a fine spacing 
in stellar parameters and interpolation errors are mod- 
est. R07 estimate that interpolation within the grid is 
consistent with non-LTE corrections computed for each 
star individually within 0.02 dex, but in most cases the 
agreement is better than within 0.01 dex. 

One important limitation of the non-LTE calculations 
employed in this work is the neglect of inelastic colli- 



sions with neutral hydrogen in the statistical equilib- 
rium calculations. They tend to bring the predicted line 
strengths closer to their LTE values, thus reducing the 
size of the non-LTE corrections. Collisions with neutral 
hydrogen can be taken into account in the non-LTE cal- 
cu lations using th e modified classical Thomson formula 
by iDrawinl (11968T) for the c r oss-se ction s. As poin t ed ou t 
bv iSteenbock fc Holwegerl (|1984n and iLambertJ (|1993f l. 
among others, Drawin's formula likely provides only an 
order of magnitude estimate, and it is therefore common 
to scale its prediction with an empirical Sh factor. It is 
expe cted that this number depends on th e type of transi- 
tion (|Lambertlll99l iBarklem et al.ll2011l ). but it is often 



12 



RAMIREZ ET AL. 



assumed constant. 

Center-to-limb variation (CLV) observations of the 
triplet can be used to constrain S h, an exercise that is 
obvio usly limited to the Sun (e.g.. lAllende Prieto et al.1 
2004). Using a three-dimensional hydrodynam i c sim - 
ulation of the solar photosphere, IPereira et al.l (2009) 
find the best agreement between model and observation 
for the CLV of the triplet if an Sh — 1 is adopted. 
iRamrrez et al.l (|2006f ). on the other hand, find that the 
three lines of the triplet in the spectrum of the standard 
SDSS star BD +174708 ([Fe/H] ~ -1.7) are most con- 
sistently fitted using Sh — 10, although that work uses 
one-dimensional static models and the oxygen abundance 
inferred for this well-known moderately metal-poor star 
with that choice of Sh appears higher than usual values 
for halo stars. 

Further investigation on the impact of inelastic col- 
lisions with neutral hydrogen on non-LTE calculations 
is clearly needed to better constrain Sh as well as its 
possibly important variation across stellar parameter 
space. Note, however, that the detailed qu antum me- 
chanical calculations by IBarklem et al.l (|2011[ ) show that 
the physics of excitation by H collisions is not well repro- 
duced by Drawin's formula, which could result in poorly 
predicted relative rates for collisional transitions. 

In this work, we did not take into account collisions 
with neutral H. Therefore, there is no Sh factor associ- 
ated with our non-LTE calculations. Nevertheless, as ex- 
plained in R07, by analyzing the non-LTE oxygen abun- 
dances inferred for each of the triplet lines, we found 
that in order to obtain consistent abundances from the 
three lines (on the average for all stars) small corrections 
were needed. The non-LTE oxygen abundances from the 
7771.9 A line need to be reduced by 0.036dex while those 
from the 7775.4 A need to be reduced by 0.018 dex, keep- 
ing the abundance from the weakest line constant. The 
latter is expected to be less affected by collisions with 
neutral H, while the direction of these empirical correc- 
tions is consistent with them being due to the neglect 
of this effect, which should be more important for the 
strongest line which forms in higher atmospheric lay- 
ers with lower H densities. Thus, the impact of colli- 
sions with neutral H is reduced, although certainly not 
fully removed, with these empirical corrections, which we 
also applied in this work. Interest ingly, the more recent , 
accurate non-LTE calculations bv iFabbian et al.l ((2009) , 
which include collisions with neutral H adopting Sh = 1 , 
result in a non-LTE abundance correction for the Sun 
that is very similar to that from our work. 

We derived LTE oxygen abundances using the abf ind 
driver in MOOG. Given the strong sensitivity of the non- 
LTE effects described above to the stellar atmospheric 
parameters, and the wide range of T c g, log g, [Fe/H] val- 
ues covered by our sample, we could not rely only on a 
solar differential analysis, as it is sometimes done. Fig- 
ure [5] shows the difference in the [O/H] values derived 
using non-LTE and LTE oxygen abundances. Although 
the differences are close to zero around the solar parame- 
ters, even within a relatively small range of 100 K in T g 
and 0.1 dex in log g, the A [O/H] differences vary between 
about —0.1 and +0.1 dex. These corrections are sizable: 
non-LTE effects must be properly taking into account be- 
fore using oxygen abundances inferred from the 777 nm 
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Figure 8. Non-LTE corrections applied to the LTE [O/H] values 
inferred from our data as a function of stellar atmospheric param- 
eters. 



triplet to interpret Galactic chemical evolution patterns. 

Our LTE and non-LTE oxygen abundances, both ab- 
solute and differential with respect to the solar oxygen 
abundance, as inferred from the Oi 777 nm triplet lines, 
are listed in Table El The [O /H] errors given there cor- 
respond to only the line-to-line scatter for the three lines 
of the 777 nm triplet. The error in [O/H] due to un- 
certainties in stellar parameters is dominated by AT c s, 
which has a median value of 44 K. Propagating this error 
into our calculation of [O/H] for a Sun-like star results in 
an error of 0.035 dex. Similarly, the propagation of the 
median error in log g and [Fe/H] into the [O/H] determi- 
nation of a Sun-like star results in errors of only 0.001 
and 0.002 dex, respectively. Conservatively, hereafter we 
adopt a median [O/H] error of A [O/H] = 0.04dex due 
to uncertainties in the stellar parameters. 

Since the median line-to-line scatter for the non-LTE 
corrected [O/H] abundances is 0.02 dex, we conclude that 
our [O/H] values have a median error of about 0.045 dex. 
Estimating the error in [O/Fe] is not straightforward, 
because both the oxygen and iron abundances depend 
on the stellar parameters. Adding the oxygen and iron 
abundance errors in quadrature we obtain A [O/Fe] ~ 
0.07 dex, but this is an overestimated error. Based on 
a number of test calculations made taking into account 
the changes in the derived abundances due to artificial 
shifts in stellar parameters which correspond to their me- 
dian errors, we find that for most of our sample stars, 
A[0/Fe] ~ 0.05 - 0.06 dex. Note that we adopt [Fe/H] 
values inferred from the Fe n line analysis. For a Sun- 
like star, a change of ±44 K in T c s implies a change in 
[Fe/H] of ±0.03 dex if using Fei lines, but +0.015 dex if 
using Fen. The latter goes in the same direction with 
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Table 6 

Oxygen Abundance 



Star 


7779 A 
1 II £ 


77 7 A A 


777^ A 
II 1 O rt 


7779 A 
1 1 1 z ±\ 


■i^yj ,1N Li 1 Hi 

77 7 A A 


777^ A 
1 1 1 O 




[0/H]nlte 


HD 32071 


8.71 


8.69 


8.71 


8.47 


8.45 


8.50 


-0.08 ±0.01 


-0.18 ±0.02 


HD 59490 


8.67 


8.66 


8.71 


8.54 


8.55 


8.59 


-0.10 ±0.03 


-0.09 ±0.02 


HD 67163 


8.88 


8.87 


8.89 


8.70 


8.69 


8.73 


0.10 ±0.02 


0.05 ± 0.02 


HD 82960 


9.12 


9.09 


9.09 


8.75 


8.73 


8.76 


0.32 ±0.01 


0.09 ±0.01 


HD 130047 


8.96 


8.94 


8.94 


8.80 


8.80 


8.81 


0.16 ±0.01 


0.15 ± 0.01 


HD 144070 


9.31 


9.19 


9.12 


8.94 


8.85 


8.81 


0.42 ± 0.07 


0.21 ± 0.05 


HD 170058 


9.03 


9.04 


9.05 


8.83 


8.83 


8.87 


0.26 ±0.02 


0.19 ±0.02 


HD 171029 


8.90 


8.90 


8.96 


8.63 


8.62 


8.70 


0.14 ±0.04 


-0.00 ± 0.04 


HD 183490 


9.10 


9.09 


9.03 


8.81 


8.81 


8.80 


0.29 ±0.02 


0.15 ±0.01 


HD 213746 


8.84 


8.87 


8.88 


8.64 


8.67 


8.71 


0.08 ± 0.03 


0.02 ± 0.03 


HD 223723 


8.94 


8.90 


8.94 


8.72 


8.70 


8.75 


0.14 ±0.02 


0.07 ± 0.02 


HIP 171 


8.54 


8.44 


8.52 


8.36 


8.27 


8.36 


-0.26 ±0.02 


-0.30 ±0.02 


HIP 348 


8.73 


8.72 


8.68 


8.57 


8.56 


8.54 


-0.07 ±0.01 


-0.10 ±0.01 


HIP 394 


8.71 


8.64 


8.64 


8.42 


8.37 


8.40 


-0.11 ±0.02 


-0.25 ±0.01 


HIP 475 


8.55 


8.53 


8.53 


8.36 


8.36 


8.38 


-0.25 ±0.01 


-0.30 ±0.01 


HIP 493 


8.69 


8.67 


8.58 


8.50 


8.49 


8.43 


-0.14 ±0.04 


-0.18 ±0.03 


HIP 522 


9.11 


9.02 


8.99 


8.81 


8.75 


8.76 


0.26 ±0.04 


0.12 ±0.02 


HIP 530 


9.10 


9.06 


9.03 


8.80 


8.78 


8.77 


0.28 ± 0.02 


0.13 ±0.01 


HIP 544 


8.98 


8.93 


8.89 


8.88 


8.83 


8.80 


0.15 ±0.02 


0.18 ±0.03 


HIP 656 


8.98 


8.93 


8.87 


8.64 


8.62 


8.59 


0.14 ±0.03 


-0.04 ±0.02 



1 Tablc[6]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. 



respect to the change in [O/H], implying that the T e g 
error compensates somewhat in Fe and O separately to 
make the [O/Fe] uncertainty smaller. 

For the Sun, we derive an oxygen abundance A<j = 
8.77, with a line-to-line scatter of only 0.011 dex. After 
applying the non-LTE corrections, this value reduces to 
Aq = 8.64. Our non-LTE solar oxygen abundance is 
in good agreemen t (with in 0.05 dex) with that measured 
bv lAsplund et all (|2009f) using a variety of spectroscopic 
indicators in addition to the 777 nm Oi triplet, including 
a number of [O i] lines which are expected to be formed 
in LTE. Moreover, our solar oxygen abundance is within 

0. 1 dex of the abundan ce inferred from the triplet lines 
bv lAsplund elaTl (pOOl . 

4. OXYGEN ABUNDANCE PATTERNS 
4.1. Fine-tuning the Sample 

Even though significant improvements to the determi- 
nation of stellar parameters and oxygen abundances have 
been made since our previous publication on this topic 
(R07), inspection of abundance trends, in particular the 
[O/Fe] versus [Fe/H] relation, revealed a number of out- 
liers and small subsamples of stars that added significant 
scatter to the mean relations observed. This suggests 
that there are systematic errors beyond our control that 
are still affecting our results, albeit only for a limited 
number of stars. 

In Figure we show the [O/Fe] versus [Fe/H] rela- 
tion derived using the [Fe/H] values inferred from the 
Fe I lines only. In Figure we show the same relation, 
but for Fen. Note that Figure [9] does not discriminate 
our sample stars according to their kinematics. In both 
Figures [9^ and b the well-kno wn general trend of th e 
Galactic disk is observed (e.g.. lEdvardsson et al.l fl993). 

1. e., increasing [O/Fe] abundance ratios with decreasing 
[Fe/H], reaching a nearly constant value of [O/Fe] ~ 0.45 
at [Fe/H] < —1, where it connects with the halo. How- 
ever, the star-to-star scatter is smaller when using [Fe/H] 



values from the Fe n lines alone. In the majority of our 
sample stars, Fen is the dominant species of iron, which 
makes the Fe II lines less sensitive to uncertainties in the 
stellar parameters and possibly also to modeling errors. 
Thus, although based on a fewer number of features, 
therefore implying larger random errors, the [Fe/H] val- 
ues from the Fe n analysis are more robust and they are 
adopted hereafter as the stellar metallicities. 

Severe problems have been reported for both the iron 
and O I 777 nm triplet oxygen abundance determinations 
in cool K-type dwarf stars, particularly at high metal- 
licities. Chemical abundance measurements of stars in 
young open clusters show that for T c g < 5000 K, the iron 
and oxygen abundanc es diverge towards very high values 
as T e ff decreases (e.g..lYong et al.ll2004l ; iMorel fe Micelal 
2004; iSchuler et al.ll2006bfl . However, the precise way m 
which these abundances diverge from the expected val- 
ues seems to depend on propertie s of th e cluster such 
as metallicity and age. iRami'rej (|2008| ) showed that 
this discrepancy is unlikely to be related to uncertain- 
ties in the modeling of stellar atmospheres due to the 
simplicity in the treatment of surface convection. The 
impact of chromospheric activity, starspots, and non- 
LTE effects remains to be fully explored and it could 
help understanding these problems. In Figure Eh, stars 
with T ef f < 5100 K, \ogg > 4.4, and [Fe/H] > -0.1 are 
highlighted with large crosses (11 stars). With few ex- 
ceptions, these objects are well above the mean trend. 
Given their very uncertain atmospheric parameters and 
iron/oxygen abundances, we exclude these stars from fur- 
ther analysis. 

As noted before, our sample consists mainly of dwarf 
and subgiant stars. Nevertheless, a non-negligible num- 
ber of giant stars are also included (35, defined here by: 
T off < 5500 K and log 3 < 4.0). They are highlighted 
with large crosses in Figure EH (red in the color ver- 
sion). As will be shown later, Galactic disk stars seem 
to follow distinct abundance patterns depending on their 
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-0.2 7 b) Fe II 




-0.2 r c) exclude very cool dwarfs 




-0.2 r d) exclude giants and disk outliers 
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-0.2 r e) final ~. 
-0.4 F ................................ : 
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[Fe/H] 

Figure 9. a) [O/Fc] versus [Fe/H] for all stars in our sample; 
[Fe/H] is based on Fel lines only, b) As in a) for [Fe/H] based on 
Fell lines, c) As in b); very cool dwarf stars are marked with large 
crosses and excluded from further analysis, d) As in c); giant stars 
and objects which clearly depart from the mean thin/thick disk 
trends are marked with larger crosses and excluded from further 
analysis, e) Final [O/Fe] versus [Fe/H] data adopted in this work. 

kinematics. Detailed inspection of the giant star data in 
Figure [9ji, separating the stars into thin- and thick-disk 
members as described in Section I2.4[ shows a pattern 
nearly identical to that of the dwarf and subgiant stars, 
but shifted slightly in both [O/Fe] and [Fe/H]. Thus, 
adding the giant stars to the dwarf/subgiant data in- 



creases the star-to-star scatter in the oxygen abundance 
trends. One could attempt to put the giant star data into 
the dwarf/subgiant scale by applying constant offsets in 
[O/Fe] and [Fe/H], but this type of empirical procedure is 
risky; systematic errors do not always exhi bit such linear 
behavior. Indeed, lAlves-Brito et all (|2010D have demon- 
strated that the practice of combining dwarf/subgiant el- 
emental abundance data with giant star data can lead to 
spurious results. Given that our methods are optimized 
for the study of dwarf and subgiant stars, hereafter we ex- 
clude from our work objects on the red giant branch. We 
should note, however, that for the particular case of oxy- 
gen, the dwarf /giant discrepancies may be present when 
using the 777 nm O I triplet lines and not other features 
due to oxygen. For example, in their e xtensive abun- 
dance analy sis of stars in the local region, ILuck fe Heiterl 
(2006, 2007) find that their oxygen abundance trends for 
dwarf and giant stars are equivalent. They used only the 
[O i] line at 630.0 nm for their oxygen abundance analy- 
sis. This feature is expected to be less model dependent 
than the infrared triplet, but its analysis is complicated 
by the fact that it is a weak feature blended by a Ni line. 

Even after eliminating stars with uncertain parameters 
and abundances, or inconsistent abundance trends rela- 
tive to the bulk of our data, a few disk objects (four) 
were found as outliers relative to the main disk [O/Fe] 
versus [Fe/H] trend: HIP 4039, HIP 22060, HIP 28103, 
and HIP 49988. These objects are marked with large 
crosses in Figure IHJi (blue in the color version). It is pos- 
sible that their abundances are peculiar, but we cannot 
rule out the possibility that a cool faint companion is 
affecting their photometry while not appearing obvious 
in the composite visible spectra. Very uncertain paral- 
laxes could also be partially responsible for these outliers. 
Hereafter, these objects are also excluded from our work. 
The sample finally adopted for further analysis (consist- 
ing of 775 stars, or 94 % of all our sample stars) is shown 
in Figure The overall pattern of increasing [O /Fe] 
with decreasing [Fe/H] for disk stars now appears very 
clean. Its detailed structure is discussed below in Sec- 
tions |0] to H3] The few objects with low [O/Fe] - 0.0 
at [Fe/H] ~ — 1.0 are halo stars whose somewhat peculiar 
abundances are discussed in Section 14.61 

4.2. Oxygen Abundance of Solar-metallicity Stars 

As mentioned before, one of the most important im- 
provements made with respect to our previous work 
(R0 7) is the use of the u pdated IRFM T cS scale 
by iCasagrande et alj (|2010f) . The average +100 K 
offset in T e g between tha t T e g scale and that by 
iRamirez fe Me lcndcz (2005a;Q), which was used in R07, 
has resulted in important offsets in [O/H] and [Fe/H]. 
This is particularly interesting when looking at the 
[O/Fe] of stars with near solar metallicity ([Fe/H] = 0). 
In R07, these stars seemed to have an average [O/Fe] ~ 
0.1 (see, for example, their Figure 8), leaving the Sun as 
a somewhat peculiar star regarding its oxygen content. 
Although not completely beyond the thin-disk [O/Fe] 
dispersion at that [Fe/H], the Sun appeared to have a 
marginally low oxygen-to-iron abundance ratio relative 
to other stars with similar parameters. Using our newly 
derived iron and oxygen abundances, however, the loca- 
tion of the Sun is in much better agreement with that of 
most disk stars, as suggested by Figure [9J3. 
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The work s by IMelendez etHl (|2009t ) and 
I Ramirez et al.l (|2009f) on extremely high precision 
elemental abundances of solar twin stars, objects 
which have stellar parameters so similar to the Sun 
that systematic errors in the abundance analysis are 
minimized, show that the Sun in fact has a slightly 
higher oxygen- to-iron abundance r atio relative to its 
twins (see also IRamirez et ail [2010). With abundance 
errors of about 0.01 dex (quantified as the standard 
error for the mean abundance ratios of small samples 
of solar twins), they find t hat the mean [O/Fe] of solar 
twins is slightly sub-solar. [Mclcndcz et al.l ([2009) derive 
an a verage [O/Fe] = —0.0 33 ± 0.011 for 11 solar twins 
while IRamirez et al.1 ((20091) find [O/Fe] = -0.015±0.006 
using 22 of those objects. Restricting our sample to 
objects with T c g = 5777 ± 100 K, log 5 = 4.44 ±0.1, and 
[Fe/H] = .0 ±0.1, which is the p ractical definition of so- 
lar twin by IRamirez et all (]2010f ). we find 47 stars with a 
mean [O/Fe] = 0.02 ± 0.04, which is consistent with zero 
within our uncertainties, but also with the slightly low 
solar oxygen abundance s uggest ed by IMelendez et al.l 
(|2009f) and IRamirez etHI (|2009l) . Mo reover, in their 
photo m etric studies of s olar tw in stars, [Mclcndc z et al 



2010), IRamirez et~aTI (|2012bD and ICasagrande et al 
2012) have concluded that the ICasagrande et al.l (2010: 



T e ff scale should be corrected by about +20 K in order 
to be in perfect agreement with the highly-precise solar 
colors derived using mo del- independent methods. An 
increase of +20 K in T e g implies a decrease of about 
0.01 dcx in [O/Fc], bringing our average solar twin 
[O /Fe] abundance ratio e ven closer to the values fo und 
bv IMelendez et al.l (|2009D and IRamirez et al.l (|2009l ). 

4.3. Disk: Kinematic Abundance Trends 

As demonstrated by a number of works, oxygen abun- 
dance patterns are different for nearby disk stars with 
dissimilar kinematic properties. Galactic disk stars with 
cold (warm) kinematics are associated with the thin 
(thick) disk. It is generally accepted that the thin-disk 
stars are, as a sample, more iron-rich than thick-disk 
stars, and have lower [O /Fe] rati os relative to thick 
disk s tars of similar [Fe/H] ( e.g.. iGratton et all 
20001: IProchaska et all 120001: iTautvaisiene et al/ 



1996 



2001 



Bensbv et al.l 120041: iZhang fe Zhaol 120061: IRamirez et all 
20071) . While the latter seems to be the case below 



[Fe/H] ~ —0.2, for higher [Fe/H] values, thin- and thick- 
disk star [O/Fe] versus [Fe/H] trends appear to converge 
to a common relation. 

The difference in [O /Fe] abundance ratios of thin- and 
thick-disk stars described above has been attributed to 
their separate origin. Thick-disk stars are thought to 
have been born from material enriched mainly by Type II 
supernovae (SNII) yields , which are high in O/Fe (e.g., 
IWooslev fc Weaverlfl995h . Thin-disk stars, on the other 
hand, may have originated from gas contaminated also by 
Type la supernovae (SNIa ) yields, which have high Fe/O 
(e.g- Jlwamoto et aLlll999D . The latter would increase the 
[Fe/H] values while quickly decreasing the [O/Fe] abun- 
dance ratios. Thus, the knee connecting the thick disk 
0/Fe]-[Fe/H] trend with that of the thin disk at high 
Fe/H] co uld also be a signat ure of SNIa pollution, as 
argued by Be nsbv et al.l ([20041 ). Note, however, that the 
slightly decreasing [O/Fe] abundance ratios of thick-disk 
stars with [Fe/H] up to [Fe/H] ~ —0.3 were satisfactorily 



reproduced using only metallicity-dependent SNII yields 
in the simple chemical evolution model by R07. 

One of the most accepted scenarios for the formation 
of the thick-disk involves mergers with sate llite galaxies 
early i n the history o f the Gal actic disk (e.g.JQuinn et all 
19931: IWalker et al.l fl996t IVelazauez fc White! Tl999f 



Abadi et a l. 2003; Brook et ai.ll2004f ). These events could 



include the following non-exclusive processes: 1) stars 
from an original disk are heated by interactions with the 
satellites, which perturb their orbits and increase their 
eccentricities, 2) stars formed within the merging galax- 
ies are trapped by the potential of the Milky Way's disk 
as they came close to the plane, and 3) new stars form 
from freshly accreted gas. At these early times, no sig- 
nificant contribution to the chemical enrichment of the 
interstellar medium by SNIa occurred, leaving SNII, the 
end products of massive stars' evolution, as the main 
polluters. Once settled, the remaining gas from the orig- 
inal disk is flattened by the Galactic potential and disk 
rotation, from where the present day thin-disk stars are 
born. This could have happened a few Gyr after the for- 
mation of the thick disk. By this time, low-mass stars 
had sufficient time to evolve, become white dwarfs, and 
explode as SNIa (if in a binary system, under the right 
conditions to accrete enough mass, as the most accepted 
theory for SNIa explosions suggests). The large amounts 
of iron introduced into the interstellar medium by SNIa 
quickly decreased the [O/Fe] abundance ratios, explain- 
ing the relatively steep decline of [O/Fe] with [Fe/H] for 
thin-disk stars. 

Our data appear to support the general description 
given above (Figures [TUb and b), although the thin/thick 
disk separation is not as obvious as previously reported, 
but note that our kinematic membership criterion is not 
very strong, as we adopt P\ > 0.5 for thin-disk stars and 
P2 > 0.5 for thick-disk stars. As a sample, kinematically- 
selected thin-disk stars do have lower [O/Fe] ratios than 
thick-disk stars below [Fe/H] ~ —0.2, but a number of 
these stars have [O/Fe] abundance ratios very similar to 
those of thick-disk stars of similar metallicity. In fact, 
below [Fe/H] ~ —0.5, the number of (kinematically- 
selected) thin-disk stars with low [O/Fe] seems to be 
about the same as that of thin-disk stars with high 
[O/Fe]. A similar observation can be made regarding 
kinematically-selected thick-disk stars, i.e., there are a 
good number of those objects with [O/Fe] abundance ra- 
tios very similar to those seen in the mean thin disk. 
iReddy et all (|2006l ) referred to these stars with mixed 
kinematics and abundances as "TKTA" stars. 

Note that in Figure [TU] all of the disk stars from 
our sample with reliable abundances are plotted. In 
many previous works regarding elemental abundance dif- 
ferences between thin- and thick-disk stars, stars with 
kinematics intermediate between thin- and thick-disk 
stars are excluded from the analysis and chemical evo- 
lution interpretation. This is done primarily to prevent 
kinematically-heated thin-disk stars from "contaminat- 
ing" the thick-disk star sample. In R07, for example, 
we examined the abundance patterns using only stars 
with kinematic probabilities greater than 70 % (instead 
of 50 % as in this work) . Other works have applied an 
even stronger criterion of P1/P2 > 10 for the thin disk 
(in addition to Pi > 0.5) and P2/P1 > 10 for the thick 
disk (in addition to P2 > 0.5). The thin/thick disk oxy- 
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Figure 10. a) [O/Fe] versus [Fe/H] for kinematically-selected 
thin-disk stars. The solid line was drawn by hand to trace the 
behavior of the bulk of these data, b) As in a), but for thick-disk 
stars. The dashed line was drawn by hand to trace the relation 
followed by most of these data. The solid line from panel a) is 
overplotted. 
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Figure 11. a) [O/Fe] versus [Fe/H] for kinematically-selected 
thin-disk stars (crosses) and thick-disk stars (filled circles), b) 
As is a), but for a stronger kinematic mem bership criterion. The 
solid and dashed lines are as in Figure [TOl 



gen abundance patterns that we obtain using both our 
soft and the strong kinematic membership criterion are 
shown in Figure 1111 The number of TKTA stars is re- 
duced using the strong criterion, but they do not disap- 
pear completely. We find that about 10 % of stars with 
Pi > 0.5, i.e., kinematic thin-disk stars, have thick disk 
abundances. Conversely, nearly 20% of kinematically- 
selected thick-disk stars (P 2 > 0.5) have thin-disk abun- 
dances. If we adopt the strong kinematic criterion, these 
numbers reduce to 8 % and 14 %, respectively. 

The median error in our [O /Fc] measurements is about 
0.06 dex. Assuming a normal distribution for the [O/Fe] 
errors, we expect only about 1 % of thin-disk stars to have 
thick-disk abundances, and vice versa. In a more simple 
way, one should realize that if TKTA stars were due to 
under- or over-estimated [O/Fe] values, one would expect 
to observe a non-negligible number of thin-disk stars with 
[O/Fe] significantly lower than the mean [O/Fe] trend of 



the thin disk, and, similarly, a non-negligible number of 
thick-disk stars with [O/Fe] significantly larger than the 
mean [O /Fe] trend of the thick disk. We do not see that 
in our data, which implies that errors in our abundance 
analysis can not explain the observed number of TKTA 
stars. 

Of course, the probabilistic approach that we employ 
must be taken into account when trying to interpret the 
observed fraction of TKTA stars. We can estimate the 
expected fraction of stars showing thick-disk (or halo) 
abundances to which we would erroneously assign thin- 
disk kinematics, assuming a perfect separation in the 
[O/Fe] vs. [Fe/H] plane, as n{ J2p 1 >o.s( P2 + where 
ni is the number of stars with Pi > 0.5 and the sum ex- 
tends only to those ni stars. The fraction of stars with 
thick-disk kinematics, but thin-disk abundances can be 
estimated using a similar formula: n^ 1 Yip >0 siPi+Ps)- 
Interestingly, the expected fractions of TKTA stars are 
9.2 ± 0.2% and 19.1 ± 0.6%, in good agreement with 
the observed fractions of about 10% and 20%, respec- 
tively. The error bars reflect the fact that the member- 
ship probabilities have uncertainties associated to errors 
in the U, V, W velocities, which have been propagated 
into these calculations. 

The expected and observed fractions of TKTA stars 
are very similar, implying that the existence of TKTA 
stars is a natural consequence of the overlap in the kine- 
matic distributions of thin- and thick-disk stars. How- 
ever, it should be noted that when a strong kinematic 
criterion is employed, therefore minimizing the impact of 
the overlap, the expected fractions of TKTA stars reduce 
to 2.5 ± 0.1 % and 11.3 ± 0.7%. These values are smaller 
than the observed fractions of 8 % and 14 %, particularly 
for the case of kinematic thin-disk stars. We note that 
these fractions are not very sensitive to our particular 
choice of kinematic parameters for the thin/thick disks. 
If, for exampl e, instead of using the thick disk velocity 
dispersions bylSoubiran et~a l. (2003) we adopt those by 
iRobin et all {2003), which are up to 12kms 1 larger, 
these values remain nearly unchanged; in fact they are 
different only by 0.1 and 0.3% respectively. 

We must therefore conclude that it is not possible to 
tell from the full sample of stars whether the thin- and 
thick-disk populations are fully separable in kinematics 
and abundances, because of the overlap in the kinematic 
distributions, even though the expected and observed 
fractions of TKTA stars appear to be in good agreement. 
This is because the inconsistency between the expected 
and observed fractions of TKTA stars when a strong 
kinematic separation is made weakens the assumption of 
duality. By relaxing the condition of perfect thin/thick 
separation, other possibilities to explain the variety of 
chemo-dynamical properties of solar neighborhood stars 
can be explored. 

It is possible that stars born in the thin-disk, with 
thin-disk abundances, have been heated by secular inter- 
actions, for example via collisions with molecular clouds 
or spiral arms, which would naturally explain the exis- 
tence of kinematic thick-disk stars with thin-disk abun- 
dances. The existence of kinematic thin-disk stars with 
thick-disk abundances, however, cannot be explained by 
the same types of processes. Perhaps this secular heat- 
ing is the reason why we have about twice as many thick- 
disk stars with thin-disk abundances relative to thin-disk 
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Figure 12. [O/Fc] versus [Fe/H] relation for stars with kinematics 
intermediate between those of thi n- a nd thick-disk stars. The solid 
and dashed lines are as in Figure ITul 

stars with thick-disk abundances in our sample (relative 
to the total number of stars in each population). An- 
other mechanism may mix equally these two populations 
in kinematics and abundances. 

Admittedly, our sample has largely unknown selection 
functions. It was constructed by collecting high-quality 
data from a number of sources. Thus, we cannot guar- 
antee that the fractions of TKTA stars quoted before 
are representative of a volume-limited solar neighbor- 
hood sample. Nevertheless, since these objects appear 
to have been excluded from most previous works, those 
fractions are likely lower limits to the real frequencies of 
TKTA stars. 

4.4. Disk Stars with Intermediate Kinematics 

Stars with kinematics intermediate between those of 
thin- and thick-disk stars, defined here as non-halo stars 
(P 3 < 0.5) with Pi < 0.7 and P 2 < 0.7, appear to have 
an [Fe/H] distribution slightly more similar to that of 
the thin disk than the thick disk, i.e., they are, as a 
sample, more metal-rich than the kinematically warmer 
stars. However, these "intermediate kinematic" (IK) 
stars do not have a metallicity distribution intermediate 
between those of the thin- and thick-disk stars. More- 
over, their [O/Fe] abundance ratios do not reside in the 
intermediate region, as shown in Figure [151 In fact, these 
stars seem to cluster around the thin-disk or the thick- 
disk oxygen abundance trends, but mostly the former, 
as expected given the larger number of thin-disk stars 
in any given volume of the solar neighborhood. Only 
at [Fe/H] ~ —0.4 there appears to be a hint that stars 
with intermediate kinematics have intermediate abun- 
dances, but with much larger star-to-star scatter than 
either thin- or thick-disk stars, which is not expected if 
the scatter is dominated by observational errors alone. 
In this region, the star-to-star scatter is about 0.15 dex 
while the median [O/Fe] error that we estimate is only 
about 0.06 dex. 

Although stars with intermediate kinematics do not 
necessarily have intermediate abundances, their exclu- 
sion reduces the number of stars in the region of the 
[O/Fe]- [Fe/H] plane between the mean thin-disk and 
thick-disk trends. These objects have been excluded from 
most previous studies which employed strong kinematic 
membership criteria, leading to a more clear, albeit arti- 
ficial, separation between the thin disk and thick disk. If 
anything, Figure [T2l should be more representative, or at 
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Figure 13. a) Histogram of A [O /Fe] (difference in oxygen-to-iron 
abundance ratio relative to the line dividing the mean thin- and 
thick-disk [O/Fe] patterns at -0.7 < [Fc/H] < -0.1) for stars with 
intermediate kinematics (thick solid line). Similar distributions are 
shown for thin-disk (dotted line) and thick-disk (dashed line) stars, 
b) The distribution of intermediate kinematics stars is compared 
to that of the thin- and thick-disk stars as a single group (dot- 
dashed line), c) The distribution of intermediate kinematics stars 
is compared to that of the combination of low [O/Fe] thin-disk 
stars and high [O/Fe] thick-disk stars. 



least a less biased representation, of the local disk stars 
than possibly biased plots such as that presented in Fig- 
ure II lb . For this work, we made a great observational 
effort to include as many of these IK objects as possible, 
in an attempt to reduce the bias in the samples used to 
investigate the nature of the Galactic disk. 

A more detailed investigation of the nature of the IK 
stars can be done using Figure [13l where the distribu- 
tion of IK stars around the line which divides the aver- 
age thin- and thick-disk [O/Fe] versus [Fe/H] trends is 
shown (solid line in Figures 113a . b, and c) along with 
those of thin-disk (Pi > 0.5, dotted line in Figure H3"k) 
and thick-disk (P^ > 0.5, dashed line in Figure IT3"a ) 
stars. The samples are restricted to the metallicity range 
—0.7 < [Fc/H] < —0.1, i.e., away from the high metallic- 
ity end where the thin- and thick-disk trends overlap, and 
also away from the low metallicity end which contains few 
thin-disk stars. The thin-disk and thick-disk peaks are 
clearly identified on the low A [O/Fe] and high A [O/Fe] 
sides, respectively. Their distributions, however, have 
extended tails towards the intermediate region, a prop- 
erty that is probably due to the inclusion of a significant 
number of IK stars in each group. The tails of these dis- 
tributions towards lower A [O/Fe] for thin-disk stars and 
towards higher A [O/Fe] for thick-disk stars, on the other 
hand, are fully consistent with Gaussian tails dominated 
by the observational errors. 
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In Figure IT5b the distribution of IK stars (solid line) is 
compared to that of all disk stars in the [Fe/H] range 
mentioned above (dot-dashed line). They are nearly 
identical, suggesting that the IK sample is indeed a bet- 
ter representation of the Galactic disk as a whole, i.e., 
a sample of disk stars without strong dual kinematic 
bias selection. On the other hand, if one neglects the 
extended tails of the thin- and thick-disk distributions 
towards higher and lower A[0/Fc], respectively, the dis- 
tribution of disk stars appears clearly double peaked, as 
shown in Figure IT3b (dotted line). This analytical distri- 
bution was derived by fitting Gaussian tails to the low 
[O/Fe] thin-disk and high [O/Fe] thick-disk data sepa- 
rately, and extending those Gaussian distributions to- 
wards the intermediate region^ The IK star sample 
does not look like this heavily biased group of disk stars 
at all, i.e., their A [O/Fe] distribution is not bimodal. 
Thus, as we have argued before, IK stars do not popu- 
late the region intermediate between those populated by 
the majority of thin- and thick-disk stars separately in 
the [O/Fe] versus [Fe/H] plane, but they do contribute 
more to the frequency of disk stars in that region relative 
to the case when a highly biased sample of disk stars, in 
terms of their kinematics, is used. 

4.5. The Thin/Thick Disk Duality 

Having established that stars with "ambiguous" kine- 
matics and abundances, i.e., TKTA stars as well as stars 
with intermediate [O/Fe] abundance ratios, are not un- 
common in the solar neighborhood, one must question 
whether models predicting (or assuming) a thin/thick 
disk duality are realistic. We should begin to consider 
alternative scenarios in which these objects are a natu- 
ral consequence of the Galaxy's evolution. In this con- 
text, it is relevant to mention the importance of ra- 
dial migration for the fo rmation of galactic disks (e.g., 
ISellwood fc Binnevll27)0l . 

Stars in the solar neighborhood could have been born 
at different Galactocentric distances, each having a dif- 
ferent chemical enrichment history, and brought to the 
solar vicinity after several billion years of Galactic evo- 
lution. Present-day elemental abundance gradients with 
Galactocentric dista nce for both oxygen and iron have 



been observed 



e^ 



l2002t [jacobso n et al.l 



Deharveng et all l2000t iFriel et al.l 
2007t iHenrv et al.ll2010D . implying 



that radial migration provides a way to explain the range 
of elemental abundances seen in the solar neighborhood. 
The ch emical evolution model by ISchonrich fc Binnevl 
(|2009al lbl). for example, predicts an [0/Fe]-[Fe/H] dis- 
tribution that is remarkably similar to that seen in our 
data (compare, for example, their Figure 8 with our Fig- 
ure [Tlj) . Their models "include radial migration of stars 
and flow of gas through the disk." Our data provide 
support to the idea that this mechanism is one of the 

13 The dotted line in Figure [T3b corresponds to a hypothetical 
distribution of disk stars obtained using a sample selection func- 
tion similar to the one employed in this work, but biasing it further 
towards the extremes of the thin- and thick-disk [O/Fe] distribu- 
tions. If our sample were volume-limited, we would expect the 
peak of the "thin-disk side" of that distribution (A [O/Fe] < 0) to 
be significantly higher than that at A [O/Fe] > 0, resulting in a 
Gaussian-like distribution with a single peak at A [O/ Fe] ~ —0.1 
and a weak A [O/Fe] > tail. The dotted line in Figure [13b is not 
meant to represent a volume-limited sample, but a group of stars 
chosen as in our work. 



most important ones for shaping the chemo-dynamical 
properties of stars in the solar neighborhood. 

Admittedly, some cosmo logical simulations (e.g., 
iKobavashi fc Nakasatol 120111 ) can probably reproduce 
the observational properties of the solar neighborhood 
equally well. A very important exercise in this case is, 
however, to establish from first principles which mecha- 
nism^) is(are) the most relevant one(s) with regards to 
those properties, and whether such models can produce 
galaxies like the Milky Way naturally rather than as one 
particular example in a multitude of test cases. 

As argued before, excluding stars with kinematics in- 
termediate between thin- and thick-disk stars does not 
affect in a dramatic manner the observed trends, since 
the intermediate kinematics stars do not all define a clear 
intermediate abundance pattern, but rather a majority 
of them tend to follow those of either the thin disk or 
the thick disk. Nevertheless, this does not imply that 
the thin and thick disk abundance patterns are fully 
separable in both elemental abundances and kinematics. 
As larger samples have been assembled, and very strong 
kinematic membership criteria are not applied, the thin 
and thick disk oxygen abundance trends have began to 
become less well defined thanks to the larger number of 
TKTA and IK stars analyzed. Our nearby star data do 
not fully support the idea of a dual thin/thick disk. In- 
ter estingly, ana l yses of [a/Fe] data from SEGUE/SDSS 
by iBovv et all (j2012bl ld) suggest that the disk chemo- 
dynamical properties could be well fitted as a single pop- 
ulation, questionin g the reality of the thick disk itself 
(|Bovv et alll2012aD . 

4.6. Two Halo Populations 

Even though the number of halo stars in our sample is 
relatively small, we can still obtain important informa- 
tion from their oxygen abundance patterns. Figure 114b 
shows the [O/Fe] versus [Fe/H] relation of halo stars. 
There is significant star-to-star scatter, much larger than 
that seen in the disk stars. Although the abundance er- 
rors are larger for this sample of more distant, fainter 
stars with significantly weaker spectral lines relative to 
the disk sample, they cannot fully explain the large scat- 
ter observed here. Interestingly, we find that most of 
the scatter is due to stars with a Galactic space velocity 
V < — 200kms~ 1 , as s hown in Figures [14b and c. 

iNissen fc Schuster! il2()ll): have shown that in the 
disk/halo transition region of the [a/Fe] versus [Fe/H] 
plane, i.e., at about —1.5 < [Fe/H] < —0.5, stars with 
halo kinematics separate into well-defined samples of low- 
a and high-a stars, based on their [a/Fe] abundance 
ratios, where the a elements in their study are Mg, 
Si, Ca, and Ti. They also find that stars with thick- 
disk kinematics in the same metallicity range follow al- 
most exactly the trend defined by the high-a halo stars, 
low- and high-a halo stars have als o separate Cu, Zn, 
and B a to Fe abundance ratio trends ()Nissen fc Schuster] 
1201 1|) . but there is no s ignificant difference rega rding 
their lithium abundances (jNissen fc S chuster 2012). Age 
and kinematic information useful to investigate the na- 
ture of these h a lo sub -populations has been provided by 
ISchuster et al.l (|2012|). Oxygen abundances for the ma- 
jority of stars in the INissen fc Schuster! (120101) s tudy have 
been recently derived bv lRamfrez et all (|2012af) . who find 
that the low- and high-a halo stars also separate in the 
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Figure 14. a) [O/Fe] versus [Fe/H] for kinematically selected 
halo stars with V > —200 km s b) As in a) for halo stars with 
V < -200kms- 1 . 

[O/Fe] versus [Fe/H] plane and that thick-disk stars and 
high-a halo stars follow nearly indistinguishable oxygen 
abundance patt erns. 

As shown by INissen fe Schuster! (|2010D , an important 
fraction of the low-a halo stars are retrograde regarding 
their Galactic rotation velocities. In Figures[T4"b and c we 
separate our halo stars in two groups of V > —200 km s _1 
and V < — 200 km s -1 . This allows to separate the ma- 
jority of stars in retrograde orbits from the rest of halo 
stars. Clearly, most of the scatter seen in our [O/Fe] 
versus [Fe/H] relation comes from stars in retrograde 
orbits. Note, however, that not all stars with retro- 
grade orbits are low [O/Fe] stars, a n observation that 
is also consistent with the work by INissen fc Schuster! 
(2010). The star with low [O/Fe] ~ -0.1 in Fig- 
ure [Tib is HIP 12294, which has uncertain astrometric 
data, and therefore a V velocity with a large error bar 
(V = -154.5±52.4kms _1 ). Excluding it from that plot, 
we find that the high-a, "high-oxygen" halo stars from 
our sample follow a tight [O/Fe] versus [Fe/H] relation 
such that [O/Fe] is nearly constant at [O/Fe] ~ 0.45 from 
[Fe/H] ~ —2.0 to [Fe/H] ~ —0.7, from where it decreases 
smoothly down to [O/Fe] ~ 0.2 at [Fe/H] ~ -0.2. The 
star-to-star scatter of this relation is 0.065 dex, which is 
only slightly larger than the dispersion predicted by our 
observational uncertainties, leaving little or no room for 
cosmic scatter within this group of stars. 

Thus, our oxygen abund ance results for halo stars sup- 
port the observations by INissen fc Schuster! (|2010f ) re- 
garding a-element abundances, and therefore also sug- 
gest that there is a certain degree of heterogeneity in 
the abundances a n d ki nematics of nearby halo stars. 
Niss en fc Schuster! (l2010h argue that the halo is com- 
posed of two discrete populations, and that perhaps the 
halo stars with low a-element abundances (and there- 
fore those with low [O/Fe] abundance ratios in our 
work) could have been accreted from dwarf satellite 
galaxies, most notably the present-day globular cluster 
u) Cen, whereas the high-a stars were born within the 
Galaxy (see also INissen fc Sc huster 2011; Schuster et al.l 



l2012f ). Note, however, that the wCen connection is 
questionable considering the complexity of t hat pa rtic- 
ular globular cluster whereas [Ramirez et al.l ()2012a[ ) ar- 
gue thjt_a_neast_two ofthe low-a field halo stars from 
the INissen fc Schuster! (|2010f ) sample were likely born in 
globular clusters, a claim based on the very low oxygen 
abundances and very high sodium abundances of those 
two objects. 

A dual scenario seems to explain our nearby star halo 
data. As the halo was forming, stars from the Galaxy's 
building blocks, some of which may be the present-day 
dwarf satellites, were captured by the early Milky Way 
into highly eccentric orbits, including retrograde orbits. 
The chemical composition of these objects reflects that of 
their parents, i.e., systems with slow star formation rate 
which lead to low [O/Fe] as well as low [a/Fe] abundance 
ratios at relatively low [Fe/H] because the contribution of 
SNII is reduced. Observations of halo streams and tidal 
debris : heav i ly support the idea of early mergers (e.g., 
IHelmal200l IKlementl l2010t IMaiewski et al.N2012D . Gas 
belonging to the early halo itself, with more massive stars 
forming and exploding as SNII, was used to form the 
present-day high-a, high-oxygen halo stars. Since this 
gas was probably much more homogeneous than that of 
the mixture of dwarf satellite galaxies, this scenario pre- 
dicts a small star-to-star scatter in the abundance ratios 
of their surviving members. 

Other recent works point to a halo dichotomy similar 
to that described above. Us i ng da ta from SDSS, for ex- 
ample, iCarollo et al.1 (|2007l 120101 ) find evidence for an 
outer halo that has a net retrograde rotation, contrary 
to the inner halo, which is, in addit ion, slightly more 
metal-rich (see also iBeers et al.l 120121). Th is result has 
been questioned bv iSchonrich et al.1 (|2011l ) on the ba- 
sis of a possibly biased distance determination (but see 
IBeers et"aLll2012fl. On the other hand, using also SDSS 
data, IJofre fc Weiss! (|20ll find that main-sequence turn- 
off (MSTO) stars in the halo consist of a dominant pop- 
ulation that formed quickly, but they find evidence for 
a significant number of MSTO halo stars younger than 
the dominant population, and argue that the latter may 
be the remna nts of early accretion of external galaxies. 
Interestingly, iSchuster et al.l (|2012D find that the low-a 
halo stars are about 2-3 Gyr younger than the high-a 
members, supporting this idea. 

4.7. Chemical Tagging in the Galactic Disk 

Instead of using kinematics to separate the stars into 
thin- and thick-disk objects, as described in previous Sec- 
tions, and as illustrated by Figures [T5b and b, one could 
use the elemental abundances to identify the stars as 
members of one of these two groups. In general, this 
approach is referred to as chemical tagging. The dif- 
ference between panels a and b in Figure [15] is that 
the latter (panel b) corresponds to the strong kinematic 
membership criterion (i.e., P\/Pi > 10 for the thin disk, 
P2/P1 > 10 for the thick disk), and therefore excludes a 
significant number of stars from the analysis. 

In Figure[16j Galactic disk stars have been divided into 
two groups according to their [O /Fe] abundance ratios for 
any given [Fe/H]. The separation was defined in such a 
way that, on average, the two groups resemble the trends 
followed by the bulk of kinematically-selected thin- and 
thick-disk star samples. A Toomre diagram showing the 
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Figure 15. a) Toomre diagram of stars with Pi > 0.5 (crosses) 
and P2 > 0.5 (filled circles), b) Toomre diagram of stars with 
P1/P2 > 10 (crosses) and P2/P1 > 10 (filled circles), c) Toomre 
diagra m o f stars separated using chemical tagging, as shown in 
Figure IT6l The dashed lines represent the kinematic thin/thick 
disk boundary. 
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Figure 16. Disk stars are separated into two groups using their 
[O/Fe] abundance ratios. The boundary is somewhat arbitrary, 
but it has been set so that the mean trends of these groups roughly 
correspond to those of kinematically-selected thin- and thick-disk 
stars. Solid and dashed lines are as in Figure ITT1 

location of these two groups of stars is plotted in Fig- 
ure [TSt. 

Obviously, chemical tagging of Galactic disk stars does 
not result in a perfect kinematic separation, although 
average tendencies are detected such that low- [O/Fe] 
stars tend to have smaller \V\ and (U 2 + W 2 ) 1 ^ 2 ve- 



locities. High- [O/Fe] stars, on the other hand, rotate 
slower around the center of the Galaxy (i.e., have more 
negative V velocities) and their orbits reach higher al- 
titudes from the Galactic plane. However, we must not 
neglect the fact that there are many low- [O/Fe] stars 
with warm kinematics, as well as many kinematically 
cold high- [O/Fe] stars, reminiscent of TKTA stars. 

The dashed lines in Figure [15] roughly represent 
the boundary between kinematically-selected thin- and 
thick-disk stars. The region contained between the two 
dashed lines has stars of both populations, but the re- 
gions outside of this area are unique to members of one 
group or the other. By counting the number of stars on 
the left and right hand sides of this boundary, we can 
determine the fraction of stars with low- [O/Fe] (high- 
[O/Fe]) and warm (cold) kinematics. We find those num- 
bers to be 16% and 13%, respectively. Observational 
errors alone allow these values to be only about 2 %. 

Both the kinematic and chemical tagging methods of 
separating disk stars into two groups end up with a sig- 
nificant number of stars with ambiguous kinematics and 
abundances, a number that is difficult to explain with our 
estimates of the observational errors. A simple picture 
of kinematically-selected thin-disk stars with low [O /Fe] 
and thick-disk stars with high [O/Fe] seems to no longer 
apply. 

4.8. Trends with Age 

The age distributions of thin- and thick-disk stars, se- 
lected using kinematics, are shown in Figures [TTk and b. 
The latter corresponds to the case of strong kinematic 
separation. In general, kinematically-selected thick-disk 
stars are older than thin-disk members. However, there 
is considerable overlap between the two groups when all 
stars are plotted, and some overlap when a strong kine- 
matic separation criterion is applied. In Figure [T7] (and 
hereafter), only our sample stars with ages greater than 
their 3 a error are used. 

Kinematically-selected thin- and thick-disk stars do 
have different age distributions, with the former be- 
ing younger than the latter, but we do not find in our 
data convincing evidence for a sharp boundary between 
the two groups. When a strong kinematic criterion is 
adopted, the analysis of a relatively small sample of 
disk stars may lead one to conclude that there are no 
kinematically-selected thick-disk stars below a certain 
age. Note that there are only a handful of thick disk stars 
younger than about 6-8 Gyr in Figure 117b . which, due 
to their younger ages and the fact that they are mainly 
solar-type stars, are probably close to the main-sequence 
and have relatively large age errors. Thus, these stars 
could have been excluded in previous works because of 
that, leading to a sharp lower limit for the age distribu- 
tion of thick-disk stars, which may be, however, spurious. 

In Figure [T7b. the age distributions of stars separated 
according to their [O/Fe] abundance ratios, as in Fig- 
ure [TBI are shown. Very few stars with [O/Fe] > 0.2 and 
younger than 5 Gyr are found. Although the majority of 
high [O/Fe] stars are older than about 10 Gyr, a number 
of these stars are found with ages in the range from 5 to 
10 Gyr. Also, note the presence of a few stars with low 
[O/Fe] ~ +0.1 and very old ages. 

An overall trend of increasing [O /Fe] with older age is 
clearly observed in Figure [TT1 but once again it does not 
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Figure 17. a) [O/Fe] versus age for disk stars separated using 
their kinematics (Pi > 0.5 for the thin disk and P2 > 0.5 for 
the thick disk), b) As in panel a) but for a strong kinematic 
criterion (P1/P2 > 10 for the thin disk and P2/P1 > 10 for the 
thick disk), c) [O/Fc] versus ag e for disk stars separated using 
chemical tagging, as in Figure [Tm 

separate perfectly thin and thick disk stars, regardless of 
how these groups are defined, i.e., using their kinematics 
or with chemical tagging. There is also significant scatter 
in the [O /Fe] versus age relation, which could in part be 
due to observational errors, but the presence of a non- 
negligible number of thin-disk stars with very old ages as 
well as young kinematic thick-disk stars suggests that a 
sharp thin/thick disk separation is unrealistic. 

The velocity dispersion of disk stars increases with age, 
as shown for the W component in Figure HU This result, 
which is based on all disk stars with reliable ages, is in 
excellent agreement wi th that found by, e.g., the GCS 
l|Nordstr5m et all [2004) . The la star-to-star scatter in 
the W values of disk stars increases from about lOkms -1 
at 1 Gyr to 40kms~ 1 at 10 Gyr and older ages, as shown 
in the bottom panel of Figure [TBI The U and V velocity 
dispersions of disk stars in our sample also increase with 
age, as shown in Figure [T2] 

If the stars are separated as thin- and thick-disk mem- 
bers, one could conclude that the increase in the W 
velocity scatter is due to thick-disk stars only. In Fig- 
ure [20k we note that aW increases only slightly with 
age for kinematically-selected thin-disk stars, and much 
more rapidly, as well as at a higher level, for thick-disk 
objects. On the other hand, if the thin/thick disk mem- 
bership is assigned using the strong kinematic criterion, 
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Figure 18. Top panel: Galactic space W velocities as a function 
of age for all disk stars. The solid line with error bars represent the 
average value and star-to-star scatter (1 cr) at a given age. Bottom 
panel: Star-to-star scatter in W (i.e., aW) as a function of age for 
all disk stars. 
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Figure 19. Star-to-star scatter in U (top panel) and V (bottom 
panel) as a function of age for all disk stars. 

the aW of thick-disk stars appears to be nearly constant 
with age, at about 50kms _1 , as shown in Figure [20b. 
while the thin-disk aW values do not exceed 20kms _1 
at any given age. 

The results described above suggest that a proper kine- 
matic membership assignment must take into account 
the fact that the velocity dispersions of the two popula- 
tions depend on age. Although the velocity dispersions 
as a function of age could be determined with reasonable 
precision using large survey data, the uncertainties in 
the determination of ages for individual stars prevent us 
from adopting this approach. Moreover, for many stars 
no meaningful age can be derived, implying that they 
will be excluded from the analysis, potentially introduc- 
ing another sample bias. 

When disk stars are separated according to their 
[O/Fe] abundance ratios, the run of aW with age is 
very similar for both low [O /Fe] ( "thin-disk" ) stars and 
high [O/Fe] ("thick-disk") stars, as shown in Figure 120b. 
The aW values are less certain for the thick-disk at 
young ages due to the low number of young thick-disk 
stars used to compute the velocity dispersion. This ex- 
plains the "noisy" nature of the aW versus age relation 
of chemically-tagged thick-disk stars at young ages (the 
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Figure 20. a) The star-to-star scatter in W as a function of age 
for all disk stars (diamonds and solid line) is compared to that 
of kinematically-selected thin- (crosses and dotted line) and thick- 
disk (squares and dashed line) stars, b) As in panel a) but with 
a strong kinematic selection, c) As in panel a) but employing 
chemical tagging. 

same argument can be made to explain the peculiarly low 
aW value for the oldest kinematically-selected thin-disk 
stars in Figures |2"0"r and b). Admittedly, Figure l2"0b sug- 
gests that thick-disk stars still have slightly higher aW 
at any given age relative to thin-disk stars, but the differ- 
ences are not as dramatic as seen in the former case, i.e., 
when the stars' kinematic properties are used to define 
these two samples. That the latter happens is not at all 
surprising, and one could argue that this is a redundant 
statement, but the fact that chemically-tagged thin- and 
thick-disk stars show an increase in velocity dispersion 
with age nearly identical suggests again that they may 
be in fact members of a single stellar population. 

5. FINAL REMARKS 

We have measured oxygen abundances of a large sam- 
ple of nearby stars, mostly from the Galactic disk, but 
including also a number of halo members. Significant im- 
provements to the stellar parameter determination meth- 
ods have been made with respect to our previous work on 
this topic. Moreover, we have increased the number of 
stars in our study by almost 60%. In particular, we now 
include an important number of objects with kinemat- 
ics intermediate between those typically assigned to the 
thin disk and thick disk sub-samples. The inclusion of 
these objects allows us to investigate the elemental abun- 
dance patterns in a less biased way. Most previous works 
have completely ignored these objects in an attempt to 
avoid sample "contamination," but we argue that such 
approach may have led to spurious conclusions. 

The simple picture of an old, kinematically warm thick 
disk composed of stars with [O/Fe] abundance ratios 
that are always greater than those of thin-disk stars at 
[Fe/H] < —0.1, which are younger and kinematically cold 
as a sample, is challenged by our data. We find significant 
overlap in the chemical distributions of kinematically- 



selected thin- and thick-disk stars. Observational errors 
can explain a fraction of kinematically-selected thin-disk 
stars with thick-disk abundances, and vice- versa, but the 
fractions of these stars (which we refer to as TKTA stars) 
observed in our data are higher than expected on the ba- 
sis of random errors alone. Although the observed frac- 
tions of TKTA stars are compatible with expectations 
when a weak kinematic selection criterion is used (a star 
is considered to be part of a population when its member- 
ship probability Pj is higher than 0.5), this is no longer 
true when a strong kinematic criterion is adopted (i.e., 
P 1 /P 2 > 10, P 2 /Pi > 10). Thus, alternatives to the dual 
thin/thick disk picture need to be explored. 

We find that stars with thin-disk abundances but thick- 
disk kinematics outnumber by about a factor of two the 
stars with thick-disk abundances but thin-disk kinemat- 
ics. We attribute this observation to secular perturba- 
tions of "true" thin-disk star orbits, which, however, do 
not explain the thick-disk counterparts. We have looked 
in detail at stars with kinematics intermediate between 
those of the majority of thin- and thick-disk stars, which 
are probably more numerous (in relative numbers) in our 
work than in previous chemical abundance studies of the 
local Galactic disk. These intermediate kinematics (IK) 
disk objects do not have an [O/Fe] versus [Fe/H] rela- 
tion intermediate between those of the average thin- and 
thick-disk stars. Nevertheless, they populate the inter- 
mediate region better than any of those two groups sep- 
arately, but almost equally if a sample of kinematically- 
selected thin-disk stars is combined with a similarly cho- 
sen sample of thick-disk objects. This implies that our IK 
sample is more representative of the local Galactic disk 
as a whole than the heavily biased samples of thin- or 
thick-disk stars selected using strong kinematic criteria. 

The large size of our sample and our new observations 
of an important number of IK stars have allowed us to 
see a less biased representation of the chemical proper- 
ties of the local Galactic disk, albeit using only the very 
important element oxygen. The chemical evolution of 
the solar neighborhood appears more complex than pre- 
viously thought. With the analysis of high-quality spec- 
troscopic data for larger, volume-limited samples, or ex- 
tended samples of stars with very well known selection 
functions, we will be able to better trace the history of 
formation and evolution of the Galactic disk, the reality 
of which is, in fact, not too far ahead into the future. 
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